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The  minority  carrier  diffusion  length  is  a  critical 
parameter  in  the  development  of  next  generation 
Heterostructure  Bipolar  Transistors  (HBT)  and  highly 
efficient  solar  cells.  A  novel  technique  has  been  developed 
utilizing  direct  imaging  of  electron/hole  recombination  via 
an  optical  microscope  and  a  high  sensitivity  charge  coupled 
device  coupled  to  a  scanning  electron  microscope  to  capture 
spatial  information  about  the  transport  behavior  (diffusion 
lengths/drift  lengths)  in  luminescent  solid  state 
materials.  In  this  work,  a  numerical  model  was  developed  to 
do  a  multi-parameter  least  squares  analysis  of  transport 
images.  Results  were  applied  to  the  study  of  transport  in 
materials  at  the  forefront  of  device  technology  that  are 
affected  by  quantum  scattering  effects,  where  few  reliable 
experimental  measurements  exist.  The  technique  allows  for 
easy  localization  of  the  measurement  site,  broad 
application  to  a  range  of  materials  and  potential 
industrial  automation  to  aid  the  development  of  high  speed 
electronics  for  terahertz  devices. 
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I. 


INTRODUCTION 


A.  HISTORY 

With  the  emergence  of  the  transistor  in  1947  came  a 
revolution  in  military  affairs  (RMA)  that  has  been  evolving 
over  the  subsequent  59  years.  Today,  the  battlefield  and 
our  daily  lives  are  littered  with  electronic  devices  that 
do  everything  from  helping  us  to  see  in  the  dark  to 
enabling  global  communication  links.  The  transistors  that 
act  at  the  foundation  of  these  capabilities  are  able  to 
perform  faster  and  faster  every  "Moore  cycle".  As  the 
demand  for  faster  processing  in  smaller  electronics 
packages  has  grown,  electronics  makers  have  turned  to  a 
class  of  transistor  called  the  Heteroj unction  Bipolar 
Transistor  (HBT) . 

HBTs  are  transistors  in  which  at  least  one  of  the  two 
transistor  interfaces  is  formed  of  two  distinct  materials 
[1] .  The  primary  advantage  of  HBTs  is  their  greater 
emitter  efficiency,  defined  as  the  ratio  of  current 
injected  into  the  emitter  from  an  external  source  to  the 
leakage  current  from  the  base  to  the  emitter  under  active 
operation.  This  advantage  results  directly  from  the 

valence  band  discontinuity  at  the  emitter-base 
heteroj unction .  The  larger  barrier  to  minority  carrier 
injection  from  base  to  emitter  allows  for  a  substantial 
increase  in  the  permissible  doping  level  of  the  base  layer 
of  the  HBT,  which  reduces  sheet  resistance  and  allows  for 
thinner  base  layers  without  concern  of  emitter-collector 
leakage  in  the  cutoff  mode  of  operation.  These  advantages 
result  in  a  faster  base  transit  time,  defined  as  the  time 
required  for  the  emitter  injected  carriers  (minority 
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carriers  in  the  base)  to  diffuse  across  the  base  layer  to 
the  collector,  and  a  faster  switching  speed  for  the  HBT .  A 
schematic  and  energy  band  diagram  is  shown  in  Figure  1 
highlighting  the  valence  and  conduction  band 
discontinuities . 


(a)  (b) 

Figure  1.  (a)  Schematic  cross  section  of  an  HBT  structure, 

(b)  Energy  band  diagram  of  a  HBT  operated  under  active 

mode.  (From  Ref. [1] ) 

The  use  of  different  materials  to  provide  these 
junctions,  while  adding  complexity  to  the  design,  adds 
significant  power  amplification  benefits  and  switching 
speed  advantages  [3]  .  Figure  2  shows  a  history  of 
progression  of  highest  cutoff  frequency  (FT)  for  various 
transistors  over  time.  The  cutoff  frequency  is  defined  as 
the  frequency  where  current  gain  of  one  is  achieved.  It  is 
noted  that  the  maximum  working  frequency  of  these  devices, 
Fmax  ,  is  the  frequency  where  power  gain  becomes  unity  and  is 
below  Ft  [  4  ]  . 
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Figure  2 .  Trend  of  best  achieved  FT  for  various 
transistors.  (From  Ref [4]) 

In  order  to  achieve  these  frequencies  of  operation, 
manufacturers  increasingly  rely  upon  thinner,  more  heavily 
doped  materials  to  propagate  charge  quickly  and  efficiently 
across  the  base  to  activate  the  transistor  [1] ,  [2] ,  [5]  . 
These  increases  in  doping  level  and  the  shrinking  of 
relative  dimensions,  in  particular  of  the  base  layer,  have 
coupled  to  bring  about  new  and  interesting  regimes  that 
operate  on  the  edge  of  known  macroscopically  determined 
semiconductor  transport  properties.  In  order  to 
effectively  design  and  build  the  most  efficient  devices  in 
these  new  highly  doped  low  dimensional  regimes,  new 
techniques  that  can  extract  and  model  material  properties 
on  the  sub-micrometer  scale  will  be  necessary  [9]  . 
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B.  MILITARY  RELEVANCE 

High  speed  electronics  are  important  for  a  variety  of 
military  applications.  The  Terahertz  frequency  band, 
defined  as  the  frequency  range  between  300  GHZ  and  three 
terahertz,  is  being  explored  for  applications  in  medical 
diagnostic  imaging,  security  imaging,  and  high  bandwidth 
communications,  just  to  name  a  few.  A  quote  from  a  recent 
Defense  Advanced  Research  Projects  Agency  (DARPA)  Broad 
Agency  Announcement  (BAA)  highlights  what  might  be 
considered  the  most  critical  of  these  applications  for  our 
information  centric  battlespace. 

The  continuing  need  by  DoD  for  ultra-high 
bandwidth  communications  and  sensing  will  require 
electronics  that  operate  at  THz  frequencies. 

Given  that  advanced  microwave  satellite 
communication  systems  already  operate  near  60 
GHz,  the  instantaneous  bandwidth  required  to 
fully  monitor  the  battlespace  will  certainly 
exceed  300  GHz  early  in  the  next  century.  All 
communication  bottlenecks  must  be  removed  so  that 
surveillance  systems  can  relay  their  wideband 
measurements  to  other  locations  for  real-time 
analysis  [ 10 ]  . 

Conventional  electronic  sources  and  receivers  are 
limited  by  resistances,  capacitances,  and  transit  times, 
resulting  in  a  significant  attenuation  of  high  frequency 
power.  High  power  amplifiers  based  on  HBTs  are  beginning 
to  approach  the  realm  of  terahertz  oscillations  and  may 
provide  a  simple  semiconductor  based  solution  for  an 
integrated  coherent  terahertz  source  and  detector.  A 
greater  understanding  of  the  physics  of  electrical  carrier 
transport  in  these  highly  doped,  low  dimensional  structures 
will  aid  the  development  of  these  devices  and  provide 
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manufacturers  with  more  accurate  models  and  predictive 
design  tools. 

C.  THESIS  OVERVIEW 

In  this  thesis,  (sponsored  by  National  Science 
Foundation  DMR  0203397)  an  application  of  a  new  technique 
for  imaging  charge  transport  is  discussed  [6-8] .  A  study 
of  a  series  of  low  dimensional,  heavily  doped  AlGaAs/GaAs 
heterostructures  was  conducted  with  an  emphasis  on  the 
determination  of  the  diffusion  length  of  the  minority 

carriers  as  a  function  of  impurity  doping.  These  results 
showed  values  of  the  minority  carrier  mobility  that  can 
only  be  explained  with  the  incorporation  of  quantum 
mechanical  scattering  behavior  at  very  high  carrier 
concentrations.  This  appears  to  be  the  first  direct 
measurement  of  diffusion  lengths  and  minority  carrier 
mobilities  in  this  important  material  system  [9]  .  Chapter 
II  develops  the  mathematical  model  underpinning  the 
transport  of  minority  carriers  in  the  low  dimensional 

structures  of  interest.  Chapter  III  briefly  describes  the 

experimental  apparatus  and  the  technique  used  to  extract 
the  material  properties,  while  Chapter  IV  explores  the 

theoretical  limits  of  the  model  and  demonstrates 
experimental  evidence  of  those  limits.  In  Chapter  V  the 
experimental  evidence  showing  an  increase  in  minority 

carrier  mobility  in  heavily  doped  GaAs  «  1020(cm~3)  is 
presented,  and  the  results  are  discussed  in  the  context  of 
existing  theoretical  work. 
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II.  TRANSPORT  IMAGING  IN  THE  TWO  DIMENSIONAL  LIMIT 


A.  OVERVIEW 

SEM  charge  transport  imaging  combines  two  microscopes 
a  scanning  electron  microscope  (SEM)  to  provide  high 
resolution  charge  generation  and  an  optical  microscope  to 
image  the  transport  of  charge.  It  can  be  performed  in  any 
material  with  a  luminescent  signature  associated  with 
charge  recombination.  In  its  basic  operation,  non¬ 

equilibrium  minority  carriers  are  injected  into  the 
luminescent  semiconductor  material  by  the  SEM  and  the 
resulting  radiative  recombination  is  imaged  through  the 
optical  microscope  (OM)  .  Analysis  of  the  captured  image 
allows  quantitative,  localized  transport  measurements. 

One  application  for  this  technique  is  as  a  means  to 
perform  contact-free  measurements  of  minority  carrier 
diffusion  lengths.  This  is  a  key  parameter  for  many 
devices,  including  solar  cells,  photoconductors,  and  HBTs, 
as  discussed  in  Chapter  I.  More  conventional  techniques 
for  measuring  diffusion  lengths  are  generally  limited  by 
the  need  for  contacts  and  the  spatial  averaging  that  occurs 
in  macroscopic  electrical  measurements.  Transport  imaging 
can  determine  this  important  materials  parameter  directly 
from  a  single,  zero  bias  luminescent  spot  image, 
particularly  for  samples  in  the  2D  (thin  layer)  limit. 

B.  MODELING  LIMITATIONS  AND  ASSUMPTIONS 

In  specific  application  to  the  materials  of  interest 
from  Chapter  I  we  consider  the  case  where  a  thin  sample  is 
doped  and  the  charge  generation  rate  is  sufficiently  low  so 
that  we  are  able  to  model  the  transport  of  minority 
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carriers  in  an  approximately  constant  distribution  of 
majority  carriers.  For  example,  in  doped  samples  of  a  1  |um 
active  layer  AlGaAs/GaAs  heterostructure,  with  incident 

electrons  of  ~  15  keV,  this  means  an  electron  beam  current 
of  ~lxl0~8  A  or  less  for  material  doped  at  ~lxl018  cm-3. 

This  approximation  is  made  by  assuming  a  highest  value 
generation  rate  G  of  G  ~  Eacc/Ei,  where  Eacc  is  the  incident 
electron  energy  and  Ei  is  the  energy  required  to  produce  an 
electron/hole  pair.  For  energies  in  the  ~  5  -  40  keV 

range,  one  can  approximate  Ei  ~  3Eg  for  a  bandgap  of  Eg 
[12]  .  The  total  carrier  population  An  ~  Ap  created  then 
per  electron  is  Gx.  We  approximate  here  a  lifetime  of  x~l 
ns  and  a  probe  current  of  1  nA,  but  the  results  can  be 
scaled  accordingly.  In  this  example  our  ratio  of  resident 
majority  carriers  to  minority  carriers  is  on  the  order  of 
100.  Therefore,  our  low  injection  limit  is  valid.  The 

generation  volume  radius  for  the  electrons  was  approximated 
from  the  model  of  Kanaya-Okayama  as  ~  1.5  |mn  in  GaAs  at  30 
keV,  with  a  hemispherical  generation  volume  [13] .  For  more 
heavily  doped  materials,  or  shorter  lifetime  materials,  one 
could  use  higher  probe  currents .  Transport  imaging  can  be 
performed  outside  these  limits  with  more  sophisticated 
modeling,  but  we  will  restrict  ourselves  to  the  low 
injection  case  for  the  analysis  that  follows. 

C.  MATHEMATICAL  MODEL  DEVELOPMENT 

For  cases  where  the  diffusion  length  is  comparable  to 
or  greater  than  our  system  resolution,  diffusion  of  the 
minority  carriers  will  broaden  the  luminescent  spot.  The 
extent  of  optical  emission  then  becomes  a  function  of 
minority  carrier  diffusion  length  and  the  diffusion  length 
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can  be  directly  extracted  from  the  optical  emission  image. 
This  approach  cannot  be  easily  applied  to  bulk/thick 
samples  due  to  the  generation  volume  created  by  incident 
electrons  and  the  relatively  weak  dependence  of  the 
minority  carrier  distribution  on  diffusion  length  in  3D. 
However,  since  many  new  materials  and  devices  utilize 
primarily  thin  films,  (eg,  heterostructures,  quantum  wells 
and  specifically  the  base  regions  of  HBTs)  the  range  of 
applications  for  contact-free  diffusion  length  measurements 
is  large. 

In  order  to  extract  the  diffusion  length  from  the 
optical  image,  we  model  a  steady  state  distribution  of 
minority  carriers  created  by  a  generation  region  of  finite 
extent.  The  SEM  beam,  operated  in  a  low  injection 

configuration,  is  the  source  of  the  generation  region,  and 
our  2D  assumption  is  based  upon  the  relatively  thin  depth 
of  the  active  region  compared  with  its  extent  in  the  other 
two  dimensions. 


25  keV  e  Beam 


AIGaAs  ) 

W  25/mi 

GaAs  / 

\ 

\  Generation  Region  to m 

/  v 

AIGaAs 

.25 /mi 

e  Penetration  Radius 

GaAs  Substrate 

Figure  3.  E-Beam/Sample  interaction  schematic 


9 


The  sample  reaches  steady  state  very  quickly,  and  we 
will  describe  the  distribution  of  the  minority  carriers  in 
the  optically  active  GaAs  layer.  The  heterostructure  has 
been  modulation  doped  with  Be  (p-type) ,  and  the  minority 
carriers  are  electrons. 


Beginning  with  the  continuity  equation  for  electrons 
in  a  p-type  material: 


:d 


E  =  G  -u  +1  v.j 

,  J  n  n  n 

at  q 


where  G  is  the  generation  rate  [  y  ,  ] . 

/  cm  s 

An 

U  is  the  recombination  rate  =  —  for  low  injection. 

Tn 

dn 

—  is  the  time  rate  of  change  of  electrons  per  volume 
dt 


per  second. 


Jn  is  the  current  density  vector 


An  is  the  number  of  excess  minority  carriers  available 
for  recombination 


rBis  the  lifetime  for  electrons 

Defining  the  steady  state  current  density  for 
carriers : 

( 2 )  Jn  =  q/ujiE  +  qDnVn  [  ] 

cm  s 

is  the  mobility  of  electrons  in  GaAs  and 
E  is  the  externally  applied  electric  field. 
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D,  is  the  diffusion  coefficient  for  electrons  further 


related  to  the  Diffusion  Length  by: 
(3)  L  =  jDjn 


By  combining  the  equations  above  we  get: 


(4) 


-A+v. 


T 


n 


By  our  assumption  we  are  at  steady  state  and  therefore 
the  time  rate  of  change  of  the  electron  distribution  is 
zero.  Now,  assume  a  constant  E  field  in  the  x  direction  so 
that  E  =  qEx  ,  and  Equation  (4)  becomes: 


(4.5) 


0  =  G„ 


+  M„E 


dn  id  -2 
—  +  —  V  n 
dx  rn 


By  multiplying  through  by 


E 


and  making  the 


substitution:  S  =  /utE  ,  where  S  is  the  drift  length,  we  get 
Equation  (5) . 

(5)  W2n  +  \nx-\n  =  ^ljL 

Is  x  L1  L2 

1 .  Generation  Region 

Here  we  make  a  digression  to  discuss  the  nature  of 
Gnxn ,  the  steady  state  generation  distribution  created  by 
the  balance  between  the  continuous  SEM  injection  and 
recombination  within  the  sample.  After  Donolato  and 

Venturi  we  can  define  the  distribution  as  a  depth  dosed 
Gaussian  distribution  [15]. 


A 


g(r,z;R)  = 


go 


4  z  ^ 
kRi 


,2  o\z,R) 


R  2 no  ( z,R ) 
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Here  the  key  feature  of  the  distribution  is  the 
variance  cr  being  formed  of  two  linearly  independent 
factors : 

(7)  <j2(z,R)  =  <j2  +  a2(z,R) 

where  <j0  is  the  variance  of  the  beam  and  crs  is  the  spread 
of  the  primary  and  secondary  electrons  in  the  sample  due  to 
scattering.  z  and  R  are  the  depth  coordinate  and  the 
primary  electron  range  respectively.  R  is  a  function  of 
beam  energy  and  the  atomic  number  and  density  of  the  target 
material . 


Let  a0  be  the  measure  of  the  diameter  of  the  circle 
within  the  beam  that  contains  50%  of  the  total  beam 
current,  or  d=beam  diameter=l . 67  ao  .  From  empirical 

measurements  we  can  assume  that  the  lateral  scattering 

3 

2 

variance  a  «0.1 — [15].  R  can  be  determined  from  ref  [16] 

R 

for  a  beam  energy  of  15keV  in  GaAs  to  be  approximately 
«1.5 jum  ,  and  ~2.0 jum  for  a  25  keV  beam  energy.  Though  there 
is  variation  in  the  generation  shape  as  shown  in  Figure  3 
as  a  function  of  the  depth  (z)  we  approximate  the  variance 
as  constant  for  the  generation  region  in  our  active  layer. 
Assuming  an  average  depth  of  R/2  the  generation  variance 
becomes : 


a  (z,R)  =  ,  .36d2+0.\ 


R 


For  a  15  keV  beam  energy  and  d=\.15jum  ,  a  =  \.06jum  ,  and 
the  radius  within  which  99%  of  the  charge  will  be  generated 
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is  2yf2a  =  3jum  .  This  99%  value  will  define  the  limits  of  our 
source  region  for  numerical  integration  in  later  sections. 

We  define  our  generation  function  for  the  source  term 
and  normalize  the  output  to  1  using  Equation  (6) : 

-r'2 

(9)  g(r')  =  e1<jl 

2 .  Green' s  Function  Solution 

Returning  to  the  differential  Equation  (5)  with  the 
inclusion  of  the  source  function.  Equation  (9) : 


S  2—1 

(10)  V2n  +—n  — Tn  =  —re 
L2  L2  L2 


-r 
2  a1 


By  the  use  of  an  integration  factor  we  can  conduct  a 
change  of  variables  to  eliminate  the  nx  term  thereby  making 
Equation  (10)  into  the  Helmholtz  equation.  Substituting: 


(11)  n(x,y)  =  w(x,y)eax  into  (10)  and  combining  terms  we 

get : 


:i2 ; 


( +  w,v  +  2awx  +a2w)  +  j2  e‘U  ( +aw)~  1 


L 


v  -1 

w  =  —e 
L2 


2(7 


Combining  like  terms  and  dividing  through  by  the 
exponential : 


f  2  Sa  O 

w  +  w  +  wr 

2(2  H - — 

+ 

Cl  H - - - — 

xx  yy  x 

v  rj 

i  r  i2j 

Now  we  choose  2a  +  — =  0  or  a  = - -  in  order  to  eliminate 

L 1  2  L2 

the  derivative  term  to  obtain: 


(13) 


w„  +  w. 


yy 


S2+4L2 
4  L4 


.  Sx  -r'1 

t2  'f  2 

w=  e~L  eLry 
Is 
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which  is  nothing  more  than  the  Helmholtz  equation  where  the 


\J  +  4L 

Helmholtz  operator  is  :  V2+K 2  and  K  =  i - - - . 

2  L2 

Recognizing  that  the  Green's  Function  for  the 
Helmholtz  operator  is  the  zeroth  order  Bessell  Function  of 
the  second  kind  [16], 

(14)  G(r;r')  =  ^-Ko(k\r-r'\) 

In  v  7 

,  sJs2+4L2 

where  k  is  the  real  part  of  K:  k  = - r - . 

2  L2 

The  general  solution  to  a  Green' s  Function  problem  is 

oo  *  Sx'  —r'2 

w(x,  y )  =  J  G(r;  r')  —  e2L~  ela~  dr' 


Returning  to  the  solution  for  the  electron  distribution  ( n ) 


by  the  substitution: 


n(x,y)  =  w(x,y)eax  or  n(x,y)  =  w(x,y)e2 


Now, 


"(x’y)=^r!K° 


1  \}rl^S2+4L2, 


r  —  r'\  \e2L~  e2L*  e2a~  dr' 


By  defining  our  limits  of  integration  equal  to  the 
radius  which  contains  99%  of  the  generation  region's 
minority  carriers  as  described  by  the  Gaussian  distribution 
in  Equation  (8)  our  model  for  the  minority  carrier 
distribution  including  diffusion  and  drift  in  2D  becomes: 


1  2V2  c 

''(WV)  =  T-T  j 


S2 +41} 


S(x'-x) 

1  e  2L~  e2u~  dr' 


The  algorithm  in  Appendix  A. 5  is  a  Matlab  coded  version  of 
this  solution  using  a  double  quadrature  numerical 


integration  scheme  to  calculate  the  distribution  of 
carriers  as  they  drift/diffuse  from  the  finite  area  defined 
by  the  Gaussian  generation  function  (9)  and  bounded  by  the 

circle  of  radius  2^2a  . 


15 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


16 


III.  EXPERIMENTAL  APPARATUS 


A.  INTRODUCTION 

The  minority  carrier  distribution  in  a  semiconductor 
sample  reveals  information  regarding  the  transport 
properties  of  the  material  itself.  As  shown  in  the 
previous  chapter  the  diffusion  length  of  the  minority 
carriers  determines  the  shape  of  that  distribution.  The 
novelty  of  the  Transport  Imaging  technique  is  the 
extraction  of  the  salient  aspects  of  that  distribution  from 
an  actual  sample  in  a  controlled  and  flexible  manner.  By 
combining  the  charge  injection  and  high  resolution  electron 
imaging  capabilities  of  a  Scanning  Electron  Microscope 
(SEM)  with  the  optical  resolution  of  a  Silicon  Charge 
Coupled  Device  (CCD)  Camera,  accurate  spatial 
representations  of  these  distributions  are  captured  and 
analyzed.  A  custom  software  solution  allows  for  the 
analysis  of  these  images  and  the  fitting  of  the 
experimental  data  to  the  mathematical  model'’ s  predictions. 

B.  APPARATUS  DESCRIPTION 

The  charge  transport  imaging  instrument  combines  two 
independent  microscopes  -  a  JEOL  SEM  (See  Table  1  for 
instrument  specifications)  for  generating  charge  and  an 
optical  microscope  (OM)  for  collecting  and  imaging  the 
luminescence  emitted  from  the  recombination  process.  Using 
a  retractable  arm,  the  OM  is  placed  directly  under  the  pole 
piece  in  the  SEM,  allowing  the  electron  beam  to  pass 
through  the  center  of  the  first  optical  collecting  surface. 
The  initial  demonstration  system  modifies  an  OM  attachment 
for  the  JEOL  SEM  that  was  originally  designed  to  allow  the 
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fine  adjustment  of  sample  height  required  for  wavelength 
dispersive  X-ray  spectroscopy  (WDS) .  WDS  provides  higher 
energy  resolution  than  traditional  energy  dispersive  X-ray 
analysis  (EDX)  and  requires  sensitive  control  of  sample 
height  in  order  to  maintain  proper  conditions  for  multiple 
diffraction  angles.  OMs  designed  for  this  purpose  have 
short  focal  lengths  and  are  normally  used  in  conjunction 
with  a  lamp  source  and  low  sensitivity  near-IR/visible 
imager . 


Beam  J. 
blanker 

e  beam 


Monochromator 
and  PMT 


Optical 

Microscope 


Retractable 
CL  option 


Image 

acquisition 

and 

display 


Continuous  flow  He 


Figure  4 .  Transport  Imaging  System  Schematic 


A  schematic  of  the  system  is  shown  in  Figure  4.  In 

addition  to  the  charge  transport  imaging  microscope,  the 

instrument  is  equipped  with  standard  CL  capability  using  a 

Gatan  (formerly  Oxford  Instruments)  system  with  a  parabolic 

mirror,  k  m  monochromator  and  TE  cooled  GaAs  PMT  as  the 

detector.  Beam  blanking  capability  exists  for  future  time 

resolved  work.  Finally,  the  instrument  has  a  liquid  helium 

cooled  stage,  so  that  transport  imaging  and  conventional  CL 

can  be  performed  over  a  temperature  range  from  300  to  ~  5 

The  system  uses  continuous  flow  liquid  He,  with  the 
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K. 


sample  stage  inside  the  SEM  mounted  in  a  cold  finger 
configuration . 


Variable  accelerating  voltage;  200  to  40,000V 

Variable  probe  current;  1x10E-8  to  1x10E-12  Amps 

Maximum  sample  size  of  6"  in  any  one  dimension 

Working  distances;  8  to  48mm 

Sample  rotation;  360° 

Sample  tilting  90° 

Variable  magnification;  lOx  to  300,000x 

Maximum  resolution;  10  nm 

Secondary  and  Backscattered  Electron  detectors 

Equipped  with  EDS  capable  of  detecting  Carbon  and  forming 
X-rav  maps  of  composition;  composition  to  within  0.1  wt% 

Integrated  digital  imaging  system 

Noise  reduction  through  frame  averaging 

Image  capture  and  export  in  electronic  form  (TIFF 

Low  cost,  medium  quality  thermal  printouts 

High  quality,  medium  cost  Polaroid  type  55  film  containing 
both  negative  and  positive 

Table  1.  Table  1.  JEOL  840A  Specifications  (From  [11]) 

For  transport  imaging,  the  OM  is  used  in  a  passive 
detection  mode,  detecting  light  emitted  directly  from  the 
sample  using  a  high  sensitivity  cooled  Si  CCD  array  camera. 
The  current  camera  uses  a  2184  x  1472  pixel  array  (15  mm  x 
10  mm),  with  a  pixel  size  of  6.8  x  6.8  pm2  and  can  be  used 
for  transport  imaging  for  wavelengths  from  ~  350  to  1100 
nm.  Initial  image  processing  is  performed  using  MicroCCD,  a 
software  program  provided  with  the  CCD  camera.  Although 
further  image  and  data  processing  are  often  required  for 
individual  investigations,  we  benefit  from  excellent 
existing  image  acquisition  and  processing  capabilities, 
often  developed  to  support  astronomical  communities  using 
similar  cameras  for  low  light  imaging. 
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The  optical  microscope  insert  is  a  basic  two  lens 
system  (objective  and  eyepiece)  modified  to  allow  for 
passage  of  an  incident  electron  beam.  The  considerations, 
as  with  any  optical  microscope,  are  resolution  and 
magnification.  Estimating  the  resolution  for  incoherent 
emission  as 


(17) 


Ay 


0.61/1 

NA 


(where  Ay  is  the  spatial  resolution,  X  is  the  wavelength  and 
NA  is  the  numerical  aperture  (set  here  at  0.95  max)),  we 
find  Ay  =  0.56  |um  for  X  =  870  nm  (e.g.,  room  temperature 
emission  from  GaAs)  and  Ay  =  0.22  |um  for  X  =  350  nm  (e.g., 
emission  from  GaN)  .  The  current  magnification  of  the 

optical  system  is  ~  20x,  i.e.,  a  5  x  5  |u,rrr  area  scanned  by 

the  e  beam  creates  a  100  x  100  |jm2  area  on  the  CCD  area.  As 
mentioned,  pixel  dimensions  are  6.8  |um,  so  the  resultant 
effective  scale  for  the  final  image  is  ~0.4  |um/pixel, 
comparable  to  the  resolution  limit  for  red/near  IR  light. 
In  order  to  select  photon  emission  from  specific  regions 
within  the  sample,  appropriate  combinations  of  bandpass 
filters  are  placed  within  the  optical  path  for  wavelength 
selection.  The  filters  are  used  to  eliminate,  for  example, 
substrate  luminescence  or  to  select  the  transport  of 
interest  in  a  multilayer  sample. 

While  the  optical  resolution  limit  is  the  fundamental 
mechanical  limit  of  the  luminescence  collection,  there 
exists  a  further  analytical  bound  on  the  extraction  of 
transport  properties  related  to  the  data  extraction  method. 
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C .  DATA  EXTRACTION 

As  discussed  previously,  the  transport  property 
information  is  contained  in  the  distribution  of  minority 
carriers  at  steady  state.  This  distribution  is  directly 
linked  to  the  resulting  photon  distribution  as  captured  in 
an  image  by  our  device.  Figure  5  shows  one  such  image. 

Transport  Imaging  Luminescence  Spot 

20  ■ 

:  H 

60^| 


140  ■ 

160  H 

180  H  J 

20  40  60  80  100  120  140  160  180 

Pixel  Number  (X-axis) 

Figure  5 .  CCD  image  of  Experimental  Sample 

The  data  underlying  this  image  is  a  2x2  matrix  whose 
indices  correspond  to  pixel  number.  The  values  entered  in 
each  element  (0-10,000)  of  this  matrix  are  the  raw 
intensity  of  the  photon  emission  collected  by  the  CCD.  In 
order  to  study  the  full  extent  of  the  distribution  of 

minority  carriers  with  greatest  resolution,  we  extract  line 
segments  that  cross  the  peak  intensity  point  of  the  image. 
Though  various  methods  may  be  employed  to  select  these  data 
sets,  in  this  work  that  extraction  was  conducted  via  an 
algorithm  written  in  MATLAB  code  (See  appendix  A.l- 
"imagedatamanipulator .  m" )  .  Once  a  line  segment  is 

extracted,  it  must  be  parameterized  and  fitted  to  the  model 
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equation  developed  in  Chapter  II.  This  fitting  can  be 
accomplished  via  two  methods  with  varying  degrees  of 
flexibility  and  resolution. 

1.  Slope  Analysis  Estimation 

By  assuming  that  the  argument  of  the  Bessel  function 
is  large  compared  with  1  we  can  assume  the  distribution 
approximates  a  negative  exponential. 

—r 

With  r»Ld  :  Ko(—)  — » eLi 

so  that  the  slope  of  this  line  segment  plotted  on  a  semilog 

-1  |  e 

piot  would  be  m  = —  =  -  -  where 

Ld  ]J  inkT 

m  is  the  slope  and  r  =  *\x 1  +y2  ,  and  all  other  terms  are  as 
defined  in  Chapter  II.  Figure  6  shows  a  semilog  plot  of  a 
line  segment  extracted  from  the  image  of  Figure  5. 


Figure  6.  Semilog  plot  of  extracted  line  of  luminescence 
from  Data  Image  of  Figure  5 . 
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Figure  7  shows  the  results  of  the  slope  analysis  as 
calculated  by  the  "Slope2.m"  algorithm  of  Appendix  A. 3. 


Slope  Analysis  Estimate  of  Ld 


Figure  7.  Slope  Analysis  for  Experimental  Sample 

The  error  bars  here  are  derived  from  the  slope 
calculation  of  the  distribution  when  the  maximum  possible 
mechanical  error  limits  ( ±0Ajum  )  are  assumed  for  the  first 
and  last  points  in  the  data  sample.  This  error  analysis  is 
utilized  instead  of  a  standard  deviation  of  data  points 
from  the  linear  regression  due  to  its  physical  nature.  It 
provides  an  intuitive  link  between  the  analytical 
assessment  and  the  mechanical  limits  of  our  apparatus. 

The  benefits  of  the  slope  analysis  technique  lie  in 
its  direct  extraction  of  transport  properties  from  an  image 
with  limited  fitting  or  data  manipulation.  Depending  on 
signal  to  noise  ratio  and  sample  luminosity  it  can  provide 
knowledge  of  a  material's  diffusion  length  over  a  roughly 
3|um2  area.  This  material  property  resolution  is  limited  by 

the  optical  resolution  of  the  system  (pixel  width)  in  the  <f> 
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direction,  equal  to  OAjum ,  and  the  number  of  pixels  required 
to  conduct  the  linear  regression  analysis  in  the  r 
direction.  Figure  8  shows  a  pictorial  of  resolution 
dependence  on  data  sample  size  and  selection.  The 
resolution  listed  in  the  figure  follows  from  Figures  4-6. 


Figure  8.  Slope  Analysis  resolution 

So  in  this  case,  the  sample  resolution  of  transport 
properties  is  averaged  over  an  area  of  0.4  (im  x  6 /um  or 

2  Aiu  m  2  . 

As  mentioned,  the  sample  size  and  region  selection  is 
a  function  of  the  interplay  between  signal  to  noise  ratio, 
error  analysis  and  the  large  r  limit.  Lower  error 

estimations  require  a  larger  number  of  data  points,  while 
for  most  samples,  noise  limitations  drive  our  outer  limit 
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below  the  ideal  r/Ld>>l  limit.  In  work  done  by  M.  Talmadge 
of  Fairfield  University,  Fairfield  CT,  it  has  been  shown 
that  when  r>9Ld  the  extracted  -1/m  is  within  95%  of  the 
actual  Ld.  Figure  9  shows  the  trend  of  predicted  Ld  vs. 
actual  Ld  as  a  function  of  r/Ld[14] . 


Lexp/Lactual  VS.  x/Lactual 


x/ La  ctu  a  I 

Figure  9.  Slope  Method  Assumption  Dependence  on  large 
Bessel  Function  argument 

However,  when  Ld  is  not  known  it  is  more  difficult  to 
determine  this  confidence  factor.  Work  is  currently  being 
done  to  perfect  a  second  derivative  analysis  to  determine 
this  confidence  factor  without  a  priori  knowledge  of  the 
actual  Ld.  Additionally,  when  our  samples  are  of  low 
luminosity  and  have  short  diffusion  lengths  the  collected 
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photon  emission  does  not  possess  the  required  extent  to 
allow  data  selection  within  reasonable  limits  of  r/Ld>>l. 
In  this  case  a  drift  analysis  is  preferable  however,  it 
also  possesses  similar  limitations. 

2.  Least  Squares  2-Parameter  Fit  Analysis 

By  applying  an  iterative  least  squares  analysis  (See 

algorithm  Appendix  A.  4)  one  can  fit  the  model  prediction 
from  Equation  (16)  to  the  full  distribution  of  the 

extracted  line  of  data  and  determine  the  diffusion  length 
of  the  minority  carriers  and  the  radius  of  99%  charge 
generation.  This  analysis  technique  provides  an  additional 
understanding  of  the  generation  region  dimension  and  is 
unencumbered  by  the  limitations  of  the  r/Ld>>l  limit. 
However,  it  increases  the  area  of  the  sample  used  to 

extract  a  diffusion  length  -  effectively  reducing  the 

resolution  of  the  technique  from  «3  /urn1  to  as  large  as 
»1 6jum2  for  the  same  sample  data  from  Figure  5.  An  example 
of  a  completed  fit  is  shown  in  Figure  10,  where  n  is 

defined  as  the  radius  which  encompasses  99%  of  the  carrier 

generation  or  n  =  2^2cr  . 
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Radial  distance  from  beam  center  (cm) 


Figure  10.  Least  Squares  Model  Fit  of  Experimental  Data 

While  the  area  over  which  the  parameters  are 
determined  is  larger,  this  technique  benefits  from  the 
possibility  of  greater  accuracy  in  determination  of  those 
parameters.  Moreover,  when  the  material  studied  has  a 
diffusion  length  greater  than  our  optical  resolution,  with 
appropriately  taken  data  and  reasonable  signal  to  noise 


ratios  we  expect 

to 

extract 

diffusion  lengths 

and 

generations  region 

radii 

accurate 

to  within  ±0.1//m  . 

As 

mentioned  before. 

this 

method  is  not  limited  by 

the 

necessity  to  take 

data 

far  from 

the  generation  source. 

which  is  difficult  for  materials  of  low  luminescence. 
However,  there  are  inherent  limits  to  the  materials  and 
conditions  which  can  be  treated  with  this  analysis. 
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IV.  TRANSPORT  IMAGING  PREDICTIONS  AND  LIMITATIONS 


A.  ERROR  ESTIMATION  FOR  CURVE  FITTING  ALGORITHM 

Apart  from  noise  and  resolution  limits,  there  are  two 
primary  sources  of  error  in  this  technique;  1)  the  fit  of 
the  model  data  to  the  experimental  data,  and  2)  the 
assumptions  that  underlie  our  mathematical  model.  We  can 
treat  the  quality  of  the  model  fit  through  the  calculation 
of  a  root  mean  square  error  (RMSE) . 

The  residue  listed  in  the  legend  in  each  figure  is  the 
sum  of  the  least  squares  difference  used  to  select  the  most 
appropriate  parameter  fit.  From  this  residue  one  can 
calculate  the  RMSE  of  the  fit,  or 


RMSE  = 


vM/ 


M=1000 


where  (j) -residue-  Z(  Jpred-^ata)  and  M  is  the 

;=i 


number  of  data  points  taken  from  the  sample  and  used  in  the 
model  calculation.  This  is  a  direct  measure  of  the 
undetermined  error  of  the  fit  to  the  distribution.  Using 
this  formulation  we  routinely  achieve  very  favorable  RMSEs 

(<1(T2)  . 


As  was  noted,  our  optical  resolution  is  approximately 
OAjum ,  which  results  in  only  100  data  points  taken  over  the 
40 /urn  interval  shown  in  the  figure.  In  order  to  smooth  the 
distribution,  we  use  a  spine  interpolation  technique  to 
increase  our  data  set  by  a  factor  of  10.  For  relatively 
well  behaved  distributions,  which  these  are,  this  technique 
has  been  shown  to  not  alter  the  predictions  of  the  model, 
yet  allows  a  much  higher  confidence  in  the  curve  fit. 
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Because  the  model  distribution  is  generated  by  a 
numerical  integration  calculation  of  an  integral  expression 
and  the  parameters  are  themselves  arguments  of  non-linear 
functions,  it  is  not  readily  apparent  how  this  RMSE 
correlates  to  error  bars  in  the  parameters  themselves. 
This  analysis  can  be  done  and  is  explained  in  many  non¬ 
linear  least  squares  fitting  texts,  but  it  is  more  physical 
to  vary  diffusion  length  and  generation  radius  by  a  small 
amount  and  observe  the  resulting  magnitude  of  the  RMSE  from 


our  model  fit. 

The 

next 

series  of 

plots  show  this 

estimation . 

In  order  to  establish  a 

base  line  algorithm  precision 

we  first  produce 

a 

model 

output  for 

parameter 

values 

n  =  3.0jum;Ld  =  2.0 jum  , 

then 

allow 

the  least 

squares 

fitting 

routine  to  analyze  that  output  and  fit  it  with  the 
appropriate  parameters.  When  the  integration  step  size  is 
identical  for  the  model  produced  output  and  the  least 
squares  fitting  algorithm  the  residue  is  0  as  expected. 
This  is  shown  in  Figure  11. 
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Radial  distance  from  beam  center  (cm)  ^q-3 


Figure  11.  Baseline  calculation  error  in  Model  Fit 
algorithm  (no  variation  of  integration  step  size) 


A  slightly  more  realistic  assumption  is  that  our 
numerical  integration  step  size  is  on  the  order  of  1000 
times  smaller  than  that  of  the  experimentally  captured 
distribution.  Figure  12  shows  the  effect  of  a  difference 
in  integration  step  size  of  a  factor  of  1000  between  the 
model  produced  curve  and  the  least  squares  fitted  curve. 
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Radial  distance  from  beam  center  (cm)  x  10”3 

Figure  12.  Baseline  calculation  error  in  Model  Fit 
algorithm  (variation  of  integration  step) 

The  Residue  of  .004692  yields  an  RMSE  of  .0021  -  an 

order  of  magnitude  less  than  we  typically  see  in  our 
fitting  to  experimental  data.  We  can,  of  course,  reduce 
this  step  size  but  at  the  expense  of  integration  time  and 
given  the  very  real  effect  of  noise  in  our  collected  data, 
it  is  unrealistic  to  believe  that  we  can  achieve  greater 
precision  at  room  temperature  through  further  reductions. 

In  order  to  see  the  effect  of  parameter  errors  on  the 
Residue  we  again  employ  the  algorithm  to  produce  an  ideal 
model  distribution  and  fit  it  with  our  least  squares 
routine.  The  figures  that  follow  show  forced  errors  in  the 
model  fit  of  the  same  ideal  data  set  and  the  impact  on  the 
size  of  the  residue.  We  vary  diffusion  length  (0.1  //m  and 
0.2  /urn)  and  generation  region  radius  (0.1  jum  and0.01//m)  and 
calculate  the  RMSE  in  Table  2.  It  is  noted  that  the  RMSE 
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value  is  approximately  the  same  for  equal  variations  in 
either  parameter. 


Radial  distance  from  beam  center  (cm)  x  10^ 


Figure  13.  Error  estimation  for  OAjum  variation  in  Diffusion 

Length  RMSE=  6.2x1  (T3 


Radial  distance  from  beam  center  (cm)  x  10^ 


Figure 


14.  Error  estimation  for 
Generation  Region  radius 


OAjum  variation 

RMSE  =  5-8x10  3 


in 
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Radial  distance  from  beam  center  (cm)  x  10J 


Figure  15.  Error  estimation  for  0.01  /urn  variation  in 
Generation  Region  radius  RMSE=2.5xlO~3 

Another  reference  point  for  observed  errors  is  shown 
Figure  16. 


Radial  distance  from  beam  center  (cm)  x  10^ 


Figure  16.  Error  estimation  for  0.2 /um  variation  in 
Generation  Region  radius  RMSE=  10.5xl0~J 


in 


34 


Residue 

RMSE 

L  or  n  Variation 

.0064 

2.5x10  3 

±0.0 1  jit  in 

.0384 

6.2x10  3 

±0.1  jum 

.1110 

10.5xl0~3 

±0.2  jum 

.4763 

21.8xl0~3 

±0.5  jum 

1.6853 

41.1xl0~3 

±1.0  jum 

Table  2 .  Tabulated  Error  Estimates  for  Curve  Fits 


B.  LIMITS  OF  MODEL  ASSUMPTIONS 

A  more  fundamental  error  resides  in  the  boundaries  of 
where  our  model  assumptions  break  down,  or  where  other 
aspects  of  transport  begin  to  play  a  more  dominant  role. 
There  is  much  interesting  science  in  this  aspect  of  the 
analysis,  and  in  fact.  Chapter  V  will  focus  on  the 
interplay  of  one  such  phenomenon,  photon  recycling,  which 
at  high  doping  levels  begins  to  affect  the  luminescence 
distribution  on  a  scale  that  demands  special  treatment. 

1.  Low  Injection  Assumption 

An  important  limitation  in  our  modeling  is  our 
assumption  of  low  injection.  As  described  in  Chapter  II. B. 
for  these  samples  we  are  restricted  to  probe  currents  equal 

to  or  below  1x10  lS,4  .  Above  this  level  of  excitation  we 
significantly  alter  the  distribution  of  majority  carriers 
in  the  vicinity  of  the  generation  region  and  the 
recombination  is  no  longer  appropriately  described  as 
proportional  to  the  density  of  minority  carriers  alone.  In 
order  to  probe  this  limit  and  to  compare  the  slope  analysis 
predictions  with  the  model  fit  technique,  we  will  observe  a 
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series  of  data  taken  from  the  same  spatial  location  on  an 
AlGaAs/GaAs  heterostructure  kept  at  a  temperature  of  4.7  K, 
and  constant  beam  energy  of  25  keV  while  the  probe  current 
was  varied  from  6xlCT12 -6x10  8 ^4  .  Previous  work  within  our  lab 
has  reported  the  effect  of  increasing  SEM  probe  current  on 
the  size  of  the  luminescent  spot  [17]  .  Here  a  series  of 
images  is  presented  corroborating  this  work  and  guantifying 
the  increase  in  the  standard  deviation  of  the  generation 
distribution.  Figure  17  shows  a  schematic  of  the 
heterostructure  as  designed  by  Tom  Boone  at  Hitachi  Labs. 


Gao.6Alo.4As:  0.2pm;  p-5x1018  cm3 

electron  confinement 

Grading:  500  A 

interface  recombination 

Ga06AI04As:  0.1pm 
Na-5x1018  cm  ‘3 

PL  active  region  *870 nm 

Grading:  500  A 

interface  recombination 

Gao.6Alo.4As: 

0.2pm;  p-5x1018  cm  3 

electron  confinement 

Figure  17.  AlGaAs/GaAs  Heterostructure  design  from  Tom 
Boone  Doctoral  Dissertation  [22] 

The  study  of  these  materials  at  low  temperatures  allows  a 
greater  signal  to  noise  ratio,  which  enables  greater 
accuracy  of  the  model  fit.  Previous  work  in  our  lab  has 
shown  that  the  minority  carrier  diffusion  length  in  these 
materials  is  independent  of  sample  temperature,  and 
therefore  we  can  use  these  measurements  to  establish  a 
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baseline  of  accuracy  between  the  two  techniques  that  should 
translate  to  higher  temperatures  [23] .  Table  3  compiles  the 
salient  results  from  the  analysis  of  this  sample  at  4.7  K. 


Probe 

Current 

Ld  by  Slope 
Analysis  (|im) 

Ld  by  Model 
Fit  (|im) 

2\[2(J  Generation 
Region 

RMSE  of  Model 

Fit 

6x\0nA 

10±8 

3.9±0.5 

3.2±0.5 

25x1 03 

6x10  n  A 

2.9  ±0.2 

3.7  ±0.2 

3.8  ±0.2 

12x10  3 

6x10  10  M 

2.9  ±0.2 

4. 1  ±  0. 1 

4.2±0.1 

9x1 0^3 

6x10  ~9A 

3.4±0.2 

4.4  ±0.3 

4.5  ±0.3 

15x10  3 

6xlO~sA 

3.4±0.2 

4.7  ±0.3 

4.9  ±0.3 

15x10  3 

Table  3.  Measurement  Results  for  0.1  jum  active  layer,  Boone 


Heterostructure  #9 

2 .  Slope  Analysis  Limitations  and  the  Low  Injection 
Limit 

As  can  be  seen  from  the  tabulated  values,  increasing 
probe  current  tends  to  increase  the  effective  radius  of  the 
generation  region,  as  expected  and  reported  previously  in 
ref  [17]  .  The  diffusion  length  as  measured  by  both 
techniques  is  relatively  constant  as  a  function  of  probe 
current  below  6x10  s A  in  accordance  with  our  model 
assumptions.  Also  evident  is  the  disparity  between  the 
slope  analysis  method  and  the  model  fit.  This  difference 
is  expected  and  is  related  to  the  degree  to  which  the  slope 
analysis  limiting  assumption  of  large  Bessel  function 
argument  is  valid.  That  is,  the  slope  analysis  predictions 
assume  that  we  are  in  a  regime  where  the  Bessel  function 
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argument  — »1.  In  this  case,  measurements  were  taken  over 


a  distance  of  7-14(j,m  from  the  center  point.  Assuming  an  Ld 
actual  of  4|um  we  can  estimate  the  degree  of  disparity  that 
should  result  by  consulting  Figure  9.  Entering  the  X-axis 
with  a  value  of  7/4  or  2.5  we  extract  a  Talmadge  factor  of 
0.75,  or  we  would  expect  that  the  slope  analysis  method 
would  predict  a  value  within  75%  of  the  actual.  This 
corresponds  well  to  the  ratio  of  the  model  fit  prediction 
to  slope  analysis  prediction  3|um/4|um  -  or  0.75.  The  lowest 
probe  current  shows  the  limits  due  to  noise  in  this 
analysis.  As  expected,  the  slope  analysis  method  will  be 
impacted  more  greatly  by  poor  signal  to  noise  ratios 
because  of  its  higher  spatial  resolution. 


The  figures  that  follow  show  the  comparison  between 
slope  analysis  plots  and  model  fit  plots.  It  is 

instructive  to  observe  which  portions  of  the  model  fits 
begin  to  deviate  from  the  experimental  data  for  higher 
probe  currents.  The  trend  away  from  low  injection  can  be 
tracked  by  observing  the  deviation  in  the  "shoulder" 
regions  of  the  distributions  as  the  probe  current 
increases,  (Figures  20,21,  and  22.) 
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Figure  18.  Slope  and  Model  fit  analysis  plots  for  6x10  12 A 
probe  current  (pertinent  data  tabulated  in  Table  3) 


Figure  19.  Slope 
probe  current 


and  Model  fit  analysis  plots  for  6x10  nA 
(pertinent  data  tabulated  in  Table  3) 
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Figure  20.  Slope 
probe  current 


and  Model  fit  analysis  plots  for  6x10  10 A 
(pertinent  data  tabulated  in  Table  3) 


Figure  21.  Slope 
probe  current 


and  Model  fit  analysis  plots 
(pertinent  data  tabulated  in 


for  6x10  9 A 
Table  3) 
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Figure  22.  Slope  and  Model  fit  analysis  plots  for  6x\0  s A 
probe  current  (pertinent  data  tabulated  in  Table  3) 

From  Table  3  we  can  see  that  the  balance  between 
signal  to  noise,  and  model  limitations  place  the  best  probe 

current  for  accurate  measurement  at  6x10 ~10A  for  this  sample. 

3.  Small  Diffusion  Length  Limitations  and  the  Role  of 
the  Generation  Distribution 

A  more  substantive  measurement  limitation  for  the 
heavily  doped  materials  discussed  in  Chapter  I  is  the 
relatively  small  diffusion  lengths  that  accompany  such 
large  concentrations  of  acceptor  dopants.  While  the 
literature  predicts  an  increase  in  the  minority  electron 

mobility  in  GaAs  doped  with  Be  starting  at  5x10~I8c/iT3  ,  [23] 
the  lifetime  continues  to  trend  downward  at  a  rate  which 
overpowers  the  increase  in  mobility  and  causes  diffusion 
lengths  to  continue  to  decrease.  Just  beyond  this 

concentration  the  diffusion  length  drops  below  ljum ,  and  we 
approach  another  limit  of  our  technique.  This  limit  is 
directly  related  to  the  generation  region.  In  order  to  see 
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how  this  limit  arises,  we  again  study  our  solution  to  the 
transport  equation  for  a  Gaussian  generation  distribution. 


;  1 6 : 


n(x,y)  = 


2  jcL 


2  V2a  i 

■J*. 


V? 


+  AL- 


2L- 


\r  —  r\ 


S(x'-x)  -r’- 

|e  2l2  ela"  dr’ 


Here  the  two  terms  which  contribute  to  a  zero  E  field 
(S=0)  diffusion  are  the  Bessel  function  and  the  Gaussian 
source  function.  As  previously  discussed,  our  use  of  the 
Bessel  Equation  comes  from  well-established  differential 
equation  theory  for  solving  diffusion  equations  in  2D 
carrier  transport  and  other  disciplines  governed  by  the 
Helmholtz  Equation.  The  assumption  of  the  Gaussian 

distribution  to  represent  the  carrier  generation  region 
within  the  sample  is  based  upon  the  statistical 
interpretation  of  the  electron-electron  scattering  and  in 
limited  cases  is  backed  by  empirical  evidence 
[18] ,  [20] ,  [21]  . 

More  recently,  work  in  our  lab  has  shown  that  these 
distributions  may  not  all  fit  the  same  mathematical 
dependence  [19]  .  By  assuming  a  standard  Gaussian 

distribution,  we  may  be  neglecting  effects  of  small 
deviations  due  to  sample  geometry,  beam  inhomogeneities, 
and  possibly  other  effects  that  govern  the  sub-micrometer 
scale  granularity  of  the  minority  carrier  distribution. 
These  inaccuracies  in  our  model  will  become  more  prevalent 
when  materials  of  small  diffusion  length  are  studied. 
Moreover,  it  is  anticipated  that  when  materials  which  have 
diffusion  lengths  on  the  order  of  our  optical  resolution 
are  studied,  it  will  become  difficult  to  observe  the  effect 
of  diffusion  on  the  distribution.  At  this  limit  we  may  say 
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that  we  are  indeed  observing  the  generation  region  itself. 
As  we  approach  this  limit,  the  accuracy  of  our 
representation  of  the  interaction  region  will  become 
increasingly  important.  Any  deviation  of  the  actual 
distribution  from  our  Gaussian  model  will  be  reflected  in 
some  mixture  of  parameter  adjustments,  which  will 
unrealistically  be  portrayed  as  diffusion  length  or  sigma 
variation  by  the  fitting  algorithm. 

In  order  to  make  a  quantitative  assessment  of  this 
limit  we  will  demonstrate  a  limiting  case.  By  producing  a 
distribution  data  set  defined  as  a  pure  Gaussian  (n=3.0 jum) 
and  allowing  the  fitting  algorithm  to  fit  Equation  (16)  to 
it,  we  may  see  what  a  material  with  no  diffusion  and  a 
perfect  Gaussian  generation  distribution  might  look  like. 
Figure  23  shows  the  slope  analysis  and  model  fit  for  this 
case  and  demonstrates  that  both  methods  inaccurately 
predict  a  diffusion  length  of  0.3 jum  .  As  predicted,  the 
model  fit  compromises  the  generation  region  radius  of  the 
distribution  from  its  known  value  of  3.0 jum  in  order  to  fit 
the  data  with  a  diffusion  length  that  allows  for  the 
smallest  residue  permissible. 


Figure  23.  Pure  Gaussian  Distribution  Model  Fit 
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In  other  limiting  cases  where  the  modeled  data  set  was 
produced  with  very  small  but  non-zero  diffusion  lengths, 
the  algorithm' s  accuracy  was  directly  proportional  to  the 
step  size,  (analogous  to  pixel  size)  used  in  the  creation 
of  the  data.  This  is  consistent  with  our  prediction  that 
our  CCD  pixel  size  results  in  an  effective  lower  limit  of 
discernability  0.4|om  for  either  parameter  value.  However, 
because  the  generation  distribution  radius  does  not 
approach  this  lower  value,  it  is  effectively  only  a  limit 
for  our  determination  of  diffusion  lengths.  Compounding 
this  lower  limit  is  any  inaccuracy  in  our  assumption  of 
generation  region  form,  which  will  tend  to  distort  both  a 
and  Ld.  As  we  approach  the  regime  where  the  form  of  the 
generation  distribution  contributes  more  than  does 
diffusion  to  the  shape  of  the  overall  distribution  we 
expect  this  limitation  to  have  a  larger  and  larger  effect. 
Moreover,  the  0.4|um  error  suggested  above  is  only  accurate 
when  you  assume  a  perfect  Gaussian  generation  region.  Any 
deviation  of  the  generation  region  from  the  ideal  will  tend 
to  increase  our  baseline  error. 

In  reference  [19]  Luber  discusses  the  use  of  this  same 
Transport  Imaging  technique  as  a  means  to  more  accurately 
determine  the  interaction  region  distribution  for  materials 
of  interest.  While  a  full  quantitative  method  has  not  been 
developed,  it  is  seen  as  a  key  step  toward  the  study  of 
very  small  diffusion  length  materials  (Ld<1.0jum). 
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V.  STUDY  OF  HEAVILY  DOPED  HETEROSTRUCTURES 


A.  MOTIVATION 

As  discussed  in  Chapter  I,  faster  switching 
transistors  are  of  prime  importance  to  military 
applications.  The  devices  that  are  currently  under 
development  to  handle  this  task,  in  commercial  applications 
as  well  as  military,  are  HBTs .  The  key  device  parameter 
for  increasing  speed  and  efficiency  is  the  base  layer 
transit  time.  Many  engineering  design  approaches  are  used 
to  decrease  the  time  it  takes  electrons  to  flow  across  the 
base  layer,  but  the  efficient  use  of  these  technigues  is 
dependent  upon  the  accurate  knowledge  of  the  transport 
properties  -  minority  electron  diffusion  length,  lifetime, 
and  mobility. 


A  mathematical  description  of  the  base  transit  time 
reveals  the  importance  of  low  dimension  construction  as 
well . 


(18) 


_  W2 

=  2DP 


where  zB  is  the  base  transit  time,  W is  the  base  width,  and. 


(19) 


is  the  minority  carrier  diffusion  coefficient  as  defined  by 
the  Einstein  eguation  where  /up  is  the  minority  carrier 

mobility,  e  the  charge  on  an  electron,  k  Boltzman' s 
constant,  and  T  is  the  temperature  in  Kelvins  [1]  . 

The  appearance  of  the  base  width  as  a  squared  term 
dominates  the  trend  of  transit  time,  but  with  decreased 
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base  width  comes  the  problem  of  emitter-collector  current 
leakage.  Here,  increased  dopant  concentrations  aid  the 
reduction  of  the  base  width  by  providing  an  impediment  to 
this  current  leakage,  however,  classical  analysis  predicts 
that  increased  doping  has  a  negative  effect  on  tb  through 
its  reduction  of  / ip  and  therefore  D  [  1]  .  Classically,  one 

would  predict  a  practical  limit  to  the  concentration  of 
dopants  that  can  be  used  as  the  competing  effects  of 
reduced  Dp  and  decreased  W  interact.  However,  a  more 
detailed  analysis  reveals  a  more  complex  picture. 

B.  QUANTUM  MECHANICAL  PREDICTIONS 

An  observation  of  nonconventional  electron  current 
density  in  an  GaAs/AlGaAs  N-p-n  HBT  with  Be  base  doping  of 

6x10  lsc/«-3  led  Lyon  and  Casey  to  believe  that  some  other 
transport  mechanism  was  at  play.  They  observed  a 

collector-emitter  current  density  that  exceeded 
conventional  predictions  by  four  times  [29]  .  More 

recently,  material  growth  techniques  have  improved  and  base 
layers  are  being  produced  with  graded  doping  schemes  in  the 

low  1020  cm-3  [5]  .  At  these  levels  the  assumptions  of  the 
Boltzman  distribution  and  classical  carrier  scattering 
descriptions  may  not  be  sufficient  to  explain  carrier 
transport.  In  work  done  by  Bennett  and  Lowney  employing  a 
first  principles  quantum  mechanical  analysis  of  scattering 
mechanisms  in  heavily  doped  GaAs,  it  is  predicted  that  a 

local  minima  exists  for  electron  mobility  in  the  5jcl(T18c»r3 
regime.  This  analysis  includes  all  the  important 

scattering  mechanisms  for  the  low-field  mobilities: 
acoustic  phonon,  polar  optic  phonon,  piezoelectric,  ionized 
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impurity,  carrier-carrier,  alloy,  and  plasmon  scattering. 
The  upturn  in  the  mobility  results  from  the  dependence  of 
these  scattering  mechanisms  on  the  dopant  and  carrier 
density.  As  the  dopant  density  increases  the  average 
distance  between  holes  decreases.  This  screening  radius 
then  determines  the  upper  frequency  that  may  be  supported 
for  the  vibrational  modes  set  up  in  the  plasma  of  majority 
carrier  holes,  plasmon  cutoff  frequency  (PCF) .  As  the  PCF 
increases,  the  scattering  interaction  probability  between 
minority  carrier  electrons  and  plasmons  decreases. 
Additionally,  as  the  free  hole  concentration  increases  the 
lower  energy  bands  fill,  and  Pauli  Exclusion  Principle 
screening  becomes  important.  The  number  of  majority 
carrier/minority  carrier  scattering  events  is  reduced 
because  the  holes  are  precluded  from  changing  their  energy 
level  and  therefore  can  not  interact  [23] . 

These  results  have  been  difficult  to  reinforce 
experimentally  and  limited  direct  evidence  exists  to 
support  them  [25-28] .  A  method  to  observe  this  effect  in  a 
non-contact  manner  which  requires  little  sample  preparation 
and  is  non-invasive  would  add  to  existing  device 
diagnostics  techniques.  Transport  Imaging  provides  such  a 
solution . 

C .  EXPERIMENTAL  RESULTS 

In  order  to  observe  this  increasing  mobility  trend  a 
series  of  Be  doped  heterostructures  was  studied.  Figure  24 
shows  the  design  of  the  studied  structures. 


47 


Gao.6Alo.4As:  0.2pm;  p-5x1018  cm ~3 

electron  confinement 


Grading:  500  A  interface  recombination 

Gao.6Alo.4As:  1.0pm;  PL  active  region 

Na-3x1018  -1x102°cm'3  peak  emission  *870 nm 

Grading:  500  A  interface  recombination 

Gao.6Alo.4As:  0.2pm;  p-5x1018  cm "3 

electron  confinement 

Figure  24.  Heterostructure  design  as  grown  by  Tom  Boone, 

(From  Ref .  [ 22 ] ) 

The  Transport  Imaging  technique  was  applied  to  seven 
samples  of  differing  active  layer  doping  and  the  diffusion 
length  of  the  samples  was  extracted  using  both  the  slope 
analysis  method  and  the  two  parameter  fit.  Data  were  taken 
at  multiple  locations  on  each  sample,  on  different  days  and 
under  varying  beam  energy  and  probe  currents.  Though  some 
variation  of  parameter  value  with  location  was  noted, 
overall  the  samples  can  be  considered  to  be  very 
homogeneous,  and  these  results  to  be  representative  of  the 
average  properties.  In  order  to  test  the  predictions  of 
Dr.  Bennett,  we  required  lifetime  (x)  values  with  which  we 
could  extract  mobility  (p)  values  from  our  diffusion  length 
measurements  through  the  relationship  of  Equations  (3)  and 
(18)  .  We  coupled  independent  measurements  of  the  sample 
lifetimes  with  values  provided  by  the  grower.  Dr.  Boone.  In 
both  cases  the  measurements  were  made  by  time  resolved 
photo-luminescence  techniques.  Table  4  tabulates  the 
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initial  measured  values  for  these  samples.  The  lifetimes 
marked  with  an  asterisk  were  measured  in  2003-2004  by  Yale 
University  and  were  not  corroborated  by  our  independent  and 
more  accurate  TRPL  confirmation.  Comparisons  between  the 
Yale  reported  lifetimes  and  our  TRPL  measurements  for  the 
other  samples  showed  a  10%  overestimation  in  the  samples  of 
with  Na>xl019  cm-3. 


Sample 

Doping  [cm  3] 

x  [ps] 

Ld  [|Jm] 

cm 2 

L^J 

A2 

2.75xl018 

2050  * 

3.6±.l 

25001140 

B7 

3.75xl018 

900 

2.3±.l 

23501200 

C9 

5x1 018  I3xl018 

4800 

4.1±.l 

1333168 

D6 

3.5xl019 

95 

1.61.2 

1080012800 

E3 

4.0x1 019 

140* 

1.81.15 

920011600 

F4 

6x1 019 

116* 

1.91.1 

1240011300 

G8 

l.OxlO20  ±.lxl02° 

11 

1.71.2 

101000125300 

Table  4.  Initial  results  of  AlGaAs/GaAs  heterostructure 


study 

1 .  Initial  Observations 

The  result  for  sample  C9  is  the  only  sample  that  is 
completely  consistent  with  Dr.  Bennett's  predictions.  It 
is  also  the  sample  tested  in  Chapter  III,  and  possesses  an 
active  layer  dimension  10  times  thinner  than  the  other 
samples.  Samples  B2  and  C7  are  consistent  with  Bennett's 
trend  of  decreasing  mobility  toward  the  inflection  point  at 

a  concentration  of  5.0xl018cm~3  though  offset  by  approximately 
2 

CTYl 

1000 - .  Other  reported  data  for  mobilities  in  this  range 

V  ■  s 
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of  doping  concentrations  also  show  elevated  values  compared 
with  Bennett's  predictions  [24].  However,  the  mobilities 
for  samples  D-G  are  unrealistic,  even  if  the  trend  of 
increasing  mobility  is  evident. 

The  consistent  positive  offset  of  these  values 
suggests  a  systemic  error  or  effect  that  is  operational  in 
the  >1019crrf3  samples.  We  propose  two  reasons  for  these 
offsets  and  apply  appropriate  offsets  to  account  for  them. 


2 .  Generation  Region  Discrepancies 

As  mentioned  in  paragraph  B,  we  expect  that  as  we 
approach  the  regime  where  the  generation  region  contributes 
more  and  more  to  the  shape  of  the  extracted  curve,  we  will 
be  subject  to  limitations  due  to  inaccurate  assumptions 
about  the  generation  region.  By  studying  the  curve  fits  in 
the  higher  doped  samples  we  can  gain  some  intuition  about 
where  this  limit  may  be.  The  figures  that  follow  (Fig.  25- 
28)  are  the  best  residue  curve  fit  achieved  for  each  of  the 
samples . 


Radial  distance  from  beam  center  (cm) 


Figure  25.  Samples  A2  and  B7  best  fit  2-parameter  fit 

extractions 
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Figure  26. 


Samples  D6  and  E3  best 
extractions 


fit  2-parameter  fit 


Normalized  Experimental  vs.  Model  Fit 


Figure  27 . 


Sample  F4  best 


fit  2-parameter 


fit  extraction 


Normalized  Experimental  vs.  Model  Fit 


Radial  distance  from  beam  center  (cm)  x  10"3 


Figure  28.  Sample  G8  bes 
Assumed  reasonable  generati 

dif fus 


t  fit  (a)  algorithm  run, 
on  region  with  algorithm 
ion  length 


(b) 

fitted 
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The  generation  region  radius  has  been  shown  to 

increase  with  probe  current  and  with  beam  energy,  but  here 
we  see  a  decrease  in  the  generation  region  as  a  function  of 
dopant  concentration.  Also  noted  is  the  tendency  of  the 

model  fit  to  depart  from  the  data  set  in  the  region  of  the 

shoulders  of  the  curves  near  the  base  of  the  distribution. 

This  effectively  causes  an  overestimation  of  the  diffusion 
lengths  as  predicted  in  Section  B.  Finally,  we  can  surmise 
that  in  samples  D-G  we  are  in  the  realm  where  the  diffusion 
lengths  are  on  the  order  of  the  generation  region  radius, 
or  less,  and  approaching  a  fundamental  limit  of  our 
assumptions.  If  we  assume  that  the  produced  error  is  on 
the  order  of  that  demonstrated  with  the  pure  Gaussian  from 
Chapter  IV,  Section  C.  we  would  expect  an  overestimation  of 
the  diffusion  length  by  0.4(j,m. 

3.  Photon  Recycling  (PR) 

Another  important  and  well-documented  effect  that  must 
be  considered  is  that  of  Photon  Recycling.  The  literature 
is  replete  with  documentation  of  this  phenomenon  that 
affects  diffusion  coefficients  and  observed  lifetimes  in 
bulk  GaAs  that  begins  to  act  in  this  doping  regime  [30-32] . 
The  effect  is  treated  in  different  manners,  but 
consistently  results  in  correction  terms  being  used  to 
adjust  the  observed  diffusion  coefficient  and  total 
lifetime.  Renaud  treats  the  effect  as  an  addition  to  the 
generation  function  in  the  continuity  equation  (our 
Equation  (10))  .  He  defines  the  photon  recycling  generation 
function:  [30] 

w 

(20)  Gpr  =  K(x,x')An(x',t)dx' 
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This  represents  the  excitation  in  the  sample  with  thickness 
w  and  average  absorption  coefficient  a,  taking  into  account 
the  spectral  density  of  the  light.  xr  is  the  radiative 
lifetime  related  to  the  lifetime  we  measure  with  TRPL  by 


(21) 


where  xnr  is  the  non-radiative  lifetime.  The  real 

perturbation  comes  from  the  calculation  of  K(x,x'),  which 
is  related  as  a  series  of  exponential  integral  functions 
[30]  .  The  minority  carrier  distribution  is  then  expanded 
in  a  series  expansion  over  the  photon  recycling  source 
region  and  the  continuity  equation  is  now  adjusted  with 
each  term  possessing  a  PR  pertubation  factor  Tn,  where: 

J  w—x 

(22)  Tn= — —  [  K(x,x  +  u)undu 

n\  Tr  Jx 


Because  K  is  principally  a  function  of  exponential  integral 
functions  and  converges  quickly  to  zero  with  increasing  n, 
they  can  be  represented  by  the  spatial  average  value  (Tnj  . 
Appendix  B  lists  the  first  two  non-zero  terms  of  this 
series:  ( T0)and(T2 )  in  their  full  mathematical  form,  as  well 

as  the  exponential  integral  function.  The  resulting 
continuity  equation  is  a  modification  of  our  Equation 
(4.5)  : 


(23] 


0  =  G.  - 


1  ^ 

--ft) 


J 


j-i  di  i 

n-M„E  —  + 
dx 


f  1} 


-  +  (r2) 


V2n 


Renaud  demonstrates  good  agreement  between  his 
corrective  terms  and  behavior  of  GaAs  LEDs  and  photovoltaic 
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cells,  and  calls  for  more  study  on  small  thickness  samples. 
Badescu  states  that  while  photon  recycling  is  most  apparent 
in  bulk  samples,  there  is  a  more  pronounced  effect  in 

samples  where  the  absorption  length  La=—  exceeds  the 

a 

diffusion  length,  even  for  thinner  samples  (<l|4m)  [31]  . 

If  we  add  a  column  to  Table  4  and  populate  it  with  the 
absorption  coefficient  a  (taken  from  [33])  and  the 
absorption  length  for  each  sample  we  see  a  correlation 
between  the  departure  from  predicted  values  of  mobility  and 
the  breakpoint  where  absorption  length  exceeds  diffusion 
length . 


Sample 

Doping [ cm  3 ] 

X  [ps] 

Ld 

[  M-m] 

a 

[Om'1] 

La 

[|um] 

r  2 1 

cm 

V-S_ 

A2 

2.75xl018 

2050  * 

3.61.1 

5000 

2.0 

25001140 

B7 

3.75xl018 

900 

2.31.1 

4500 

2.2 

23501200 

C9 

5xl018  ±3x1 018 

4800 

4.1  ±  .1 

4000 

2.5 

1333168 

D6 

3.5xl019 

95 

1.6±.2 

3500 

2.85 

1080012800 

E3 

4.0xl019 

140* 

1.81.15 

3400 

2.94 

920011600 

F4 

6xl019 

116* 

1.91.1 

3000 

3.3 

1240011300 

G8 

1.0xl02°i.lxl020 

11 

1.71.2 

2300 

4.3 

101000125300 

Table  5.  HS  data  table  with  absorption  length  comparison 
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C.  CORRECTIONS  FOR  OPERATIVE  EFFECTS  AND  DISCUSSED 

LIMITATIONS 

The  photon  recycling  effect  is  dependent  upon  the 
number  of  photons  generated,  the  rate  at  which  they 
reabsorb,  but  also  on  the  rate  at  which  they  can  escape  the 
active  layer  before  creating  additional  electron-hole 
pairs.  The  first  dependencies  we  have  previously 

described,  but  now  we  must  look  at  the  index  of  refraction 
of  our  samples  and  the  corresponding  critical  angle  of 
total  internal  reflection. 


In  1976  Asbeck  reported  the  critical  angle  for 
GaAs/AlGaAs  interfaces  as  a  function  of  various  A1 
concentrations  [35] .  Interpolating  from  his  graphs  and 
confirming  with  Snell'’ s  Law,  we  arrive  at  a  critical  angle 

for  the  active  GaAs  layer  of  69.6 0  The  index  of  refraction 
varies  with  dopant  density  as  well,  but  because  the 
differences  are  small  between  AlGaAs  and  GaAs  we  can  use 
the  value  for  GaAs  as  3.59  and  for  40%  A1  concentration  in 
AlGaAs  n=3.36.  Using  these  values  to  calculate  the 
Reflectance;  [36] 


(23) 


R  = 


'  1  -  n  A " 

U  +  «J 


vi  3  36 

For  n  =  llAlGaAs  =  — —  =  .935  and  therefore  R  =  0.1%,  or  when  light 
nGaAs  3.59 

strikes  the  interface  at  an  angle  less  than  the  critical 
angle,  99.9%  will  transmit  through  to  the  AlGaAs  layer. 
Also  required  is  the  radiative  lifetime.  From  [32]  we  can 
define : 


(24) 


1 

2x1  CT10  -  Na 
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Now  to  calculate  the  correction  factors  and  apply  them  to 
our  experimental  results  for  the  samples  in  Tables  3  and  4 
we  employ  the  MATHCAD  routine  of  Appendix  B.l.  The 
detailed  calculation  sheets  are  in  Appendix  B.x  and  the 
overall  results  are  tabulated  below  in  Table  6. 


Sample 

Doping 

[cm-3] 

X  [ps] 

Ld** 

[  M-m] 

cm 2 

L^J 

tpR 

[ps] 

LdPR 

[  M-m] 

cm 2 

Her 

L v  ■  s  j 

A2 

2.75x1 018 

2050  * 

3.6±.l 

25001140 

NC 

NC 

NC 

B7 

3.75xl018 

900 

2.3±.l 

23501200 

NC 

NC 

NC 

C9 

5xl018  ±3x1 01S 

4800 

4.1±.l 

1333168 

NC 

NC 

NC 

D6 

3.5xl019 

95 

1.21.2 

600012000 

69 

1.01.2 

600012000 

E3 

4.0x1 019 

140* 

1.41.15 

560011300 

* 

OO 

1.11.15 

550011600 

F4 

6x1 019 

116* 

1.51.1 

780011000 

66  * 

1.01.1 

760011300 

G8 

l.OxlO20  f.lxlO20 

11 

1.01.2 

59000120000 

9 

0.91.2 

36000115300 

Table  6.  Tabulated  Parameters  corrected  for  Generation 


Region  error  (**)  and  Photon  Recycling  overestimation 

D.  CONCLUSIONS 

The  effect  of  the  Generation  Region  (GR)  error  plays  a 
much  stronger  role  on  the  calculated  mobility  values  than 
does  Photon  Recycling  (PR)  because  of  the  simultaneous 
effect  PR  has  on  diffusion  length  and  lifetime.  The  PR 
effect  is  seen  to  grow  as  a  function  of  doping. 

The  assumption  of  a  0.4|um  error  for  GR  is  an  estimate 
that  needs  more  refinement,  through  the  development  of  an 
analytical  assessment  of  generation  region  definition  and 
its  inclusion  in  the  numerical  integration  algorithm  of 
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normalized  electron  mobility  (/x,/V0) 


Appendix  A.  Transport  Imaging  provides  an  appropriate 
mechanism  for  this  analysis  and  should  be  pursued. 


The  mobility  values  from  the  final  corrected  column 


log  (acceptor  density  -  cm3) 

r  s  d  =  O.OIS  polyno  m  lal/transfor  m  a  +  l  o  n 

Figure  1 

IsJ 

NJ 

Figure  29.  Final  corrected  Transport  Imaging  mobility 
values  reported  (After  Bennett  [9]) 

The  local  minimum  is  clearly  demonstrated  though 

absolute  magnitude  agreement  is  not.  A  new  round  of 

experiments  is  planned  to  test  the  magnitude  relationships 
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of  the  /ur  product  through  the  measurement  of  Ldrift  by 
studying  the  distributions  as  a  result  of  an  applied  DC 
bias.  In  this  manner  a  full  distribution  fit  should  escape 
the  limitations  resulting  from  generation  region  error, 
though  the  optical  resolution  limitation  (0.4  jum)  may  still 
be  operative  in  the  samples  of  heaviest  concentrations. 

From  the  results  demonstrated  it  can  be  assumed  that 
of  the  limitations  and  constraints  inherent  in  Transport 
Imaging  the  assumption  of  a  generation  region  distribution 
has  the  largest  impact  for  measurement  of  low  diffusion 
length  materials.  It  appears  that  experimental  results  can 
be  assumed  valid  so  long  as  the  diffusion  length  measured 
is  on  the  order  of  the  generation  region  radius  (as  in 
Samples  A9-C2),  that  the  signal  to  noise  ratio  is  favorable 
(as  in  all  data  samples  shown  herein) ,  and  that  the 
diffusion  lengths  measured  are  greater  than  the  optical 
resolution  of  the  system  (0.4|um  in  these  data  samples)  . 

Several  methods  may  be  useful  to  overcome  these 
constraints  and  are  being  studied  in  our  laboratory.  They 
include  time  resolved  techniques  reminiscent  of  the  Haynes 
Shockley  experiment,  but  maintaining  the  spatial 
information  of  the  light  emission  to  great  resolution,  AC 
drift  techniques  attempting  to  generate  resonance  responses 
between  transport  properties  and  the  applied  electric 
force,  and  observation  of  the  effects  of  magnetic  fields  on 
the  flow  of  the  charge  carriers  at  the  sub  micrometer 
scale . 


58 


APPENDIX  A.  TRANSPORT  IMAGING  GRAPHIC  USER 

INTERFACE  CODE 


C:\Documents  and  Settings\FMBradley\My  Documen. . . \newf irsttry .m  Page  1 

November  16,  2005 _ 2:18:47  PM 
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function  varargout  =  newf irsttry (varargin) 

% - Log  of  program  changes - 

%29  August  -  changed  mkdir  location  for  gui  directory  creation  to  create 
%within  diffusive  data  analysis\analyzed  data. 

%Will  change  parameter  matching  to  Ld  versus  Mobility 

%19  September  -  updated  version  on  laptop  with  amended  version  from  desktop 
%12  October  -  Cleaned  up  comments  and  reorganized  order  of  functions 

%16  October  -  Reset  all  flag  values  in  each  main  analysis  functions  (af lag, gf lag,  f it 
flag)  all  tested  for  empty  variable  status  and  set  to  'O'  if  so. 

%  NEWFIRSTTRY  Application  M-file  for  newf irsttry . fig 

%  NEWFIRSTTRY,  by  itself,  creates  a  new  NEWFIRSTTRY  or  raises  the  existing 
%  singleton*. 

a 

%  H  =  NEWFIRSTTRY  returns  the  handle  to  a  new  NEWFIRSTTRY  or  the  handle  to 
%  the  existing  singleton* . 

%  NEWFIRSTTRY! 'CALLBACK' , hObject, eventData, handles, ... )  calls  the  local 
%  function  named  CALLBACK  in  NEWFIRSTTRY . M  with  the  given  input  arguments. 

%  NEWFIRSTTRY! 'Property' , 'Value' ,... )  creates  a  new  NEWFIRSTTRY  or  raises  the 

%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  newf irsttry_OpeningFunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  newfirsttry_OpeningFcn  via  varargin. 

% 

%  ‘See  GUI  Options  -  GUI  allows  only  one  instance  to  run  (singleton) . 

% 

%  See  also:  GUIDE,  GUI DATA,  GUIHANDLES 


%  Edit  the  above  text  to  modify  the  response  to  help  newfirsttry 
%  Last  Modified  by  GUIDE  v2.5  30-Mar-2005  12:09:12 


%  Begin  initialization  code  -  DO  NOT  EDIT 


gui_Singleton  =  1; 
gui_State  =  struct ( ' gui_Name ' , 

' gui_Singleton ' , 

' gui_OpeningFcn ' , 

' gui_OutputFcn ' , 

' gui_LayoutFcn ' , 

' gui_Callback ' , 
if  nargin  s  isstr (varargin(l) ) 

gui_State . gui_Callback  =  str2func (varargin) 1 }) ; 

end 


mfilename,  . . . 
gui_Singleton,  .  .  . 

@newf irsttry_OpeningFcn, 
@newfirsttry_OutputFcn, 
[]  ,  ... 

[]>; 


if  nargout 
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49  varargout {1 : nargout )  =  gui_mainf cn (gui_State,  varargin( : } ) ; 

50  else 

51  gui_mainfcn (gui_State,  varargin{ : } ) ; 

52  end 

53  %  End  initialization  code  -  DO  NOT  EDIT 

54  %Initialize  global  variables 

55  global  aflag  gflag  foldername  firsttime  fitflag; 

56  %aflag=0; 

57  %gflag=0; 

58  %fitflag=0; 

59  %firsttime=0; 

60  filetoload=[]  ; 

61 

62  % _  * 


63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 


%  -  Executes  just  before  newf irsttry  is  made  visible. 

function  newf irsttry_OpeningFcn (hObject,  eventdata,  handles,  varargin) 
%  This  function  has  no  output  args,  see  OutputFcn. 

%  hObject  handle  to  figure 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUIDATA) 

%  varargin  command  line  arguments  to  newfirsttry  (see  VARARGIN) 


%  Choose  default  command  line  output  for  newfirsttry 
handles . output  =  hObject; 


%  Update  handles  structure 
guidata (hObject,  handles); 

if  nargin  — =  3, 

initial_dir  =  pwd; 
elseif  nargin  >  4 

if  strcmpi (varargin) 1 ),' dir ' ) 
if  exist (varargin(2 ), 'dir ' ) 

initial_dir  =  varargin(2); 

else 

errordlg (' Input  argument  must  be  a  valid  directory Input  Argument  Error  if 
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95 


return 

end 

else 

errordlg (' Unrecognized  input  argument Input  Argument  Error!'); 
return; 

end 

end 

%  Populate  the  listbox 
load_listbox (initial_dir,  handles) 

%  Return  figure  handle  as  first  output  argument 


60 


C:\Documents  and  Settings\FMBradley\My  Documen. . . \newf irsttry .m 
November  16,  2005 


Page  3 
2:18:47  PM 


96  %  UIWAIT  makes  newfirsttry  wait  for  user  response  (see  UIRESUME) 

97  %  uiwait (handles .newfirsttry) ; 

98 

99  % _  4 


100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 


%  -  Outputs  from  this  function  are  returned  to  the  command  line. 

function  varargout  =  newf irsttry_OutputFcn (hobject,  eventdata,  handles) 
%  varargout  cell  array  for  returning  output  args  (see  VARARGOUT) ; 

%  hobject  handle  to  figure 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUI DATA) 


%  Get  default  command  line  output  from  handles  structure 
varargout (1)  =  handles .output ; 
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%  Callback  for  list  box  -  open  .fig  with  guide,  otherwise  use  open 


function  varargout  =  listboxl_Callback (h,  eventdata,  handles) 

%  hobject  handle  to  listboxl  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUIDATA) 

%  Hints:  contents  =  get (hobject, ' String' )  returns  listboxl  contents  as  cell  array 
%  contents {get (hobject, 'Value' ) )  returns  selected  item  from  listboxl 

get (handles . newfirsttry, ' SelectionType ' ) ; 
global  filetoload  tiftoload  ; 

if  strcmp (get (handles . newf irsttry, ' SelectionType '),' open ' ) 
index_selected  =  get (handles . listboxl, 'Value ') ; 
file_list  =  get (handles . listboxl, ' String ') ; 
filename  =  file_list{index_selected) ; 

if  handles . is_dir (handles . sorted_index (index_selected) ) 
cd  (filename) 
load_listbox (pwd, handles ) 

else 

[path, name, ext, ver]  =  fileparts (filename) ; 
switch  ext 
case  ' . fig' 

guide  (filename) 
case  '  .tif ' 

tiftoload=f ilename; 
case  ' . csv' 

filetoload=f ilename; 
otherwise 
try 

filetoload=f ilename ; 
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143  catch 

144  errordlg (lasterr , ' File  Type  Error 'modal ' ) 

145  end 

146  end 

147  end 

148  end 

149  % _  kT 


150  %  Read  the  current  directory  and  sort  the  names 

151  % - 

152  function  load_listbox (dir_path,  handles) 

153  cd  (dir_path) 

154  dir_struct  =  dir (dir_path) ; 

155  [sorted_names, sorted_index]  =  sortrows ( {dir_struct.name) ' ) ; 

156  handles . file_names  =  sorted_names ; 

157  handles . is_dir  =  [dir_struct . isdir] ; 

158  handles . sorted_index  =  [sorted_index] ; 

159  guidata (handles. newfirsttry, handles) 

160  set (handles . listboxl, ' String ' , handles . f ile_names, . . . 

161  'Value',  1) 

162  set (handles . textl, 'String' ,pwd) 

163 

164  % _  * 


165  %  -  Executes  during  object  creation,  after  setting  all  properties. 

166  function  listboxl_CreateFcn (hObject,  eventdata,  handles) 

167  %  hObject  handle  to  listboxl  (see  GCBO) 

168  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

169  %  handles  empty  -  handles  not  created  until  after  all  CreateFcns  called 

170 

171  %  Hint:  listbox  controls  usually  have  a  white  background,  change 

172  %  'usewhitebg'  to  0  to  use  default.  See  ISPC  and  COMPUTER. 

173  usewhitebg  =  1; 

174  if  usewhitebg 

175  setlhobject, 'BackgroundColor ' , 'white'); 

176  else 

177  setlhobject, 'BackgroundColor ', get (0, ' defaultUicontrolBackgroundColor ' ) ) ; 

178  end 

179 

180  % _  * 


181  %  -  Executes  during  object  creation,  after  setting  all  properties. 

182  function  edit2_CreateFcn (hObject,  eventdata,  handles) 

183  %  hObject  handle  to  edit2  (see  GCBO) 

184  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

185  %  handles  empty  -  handles  not  created  until  after  all  CreateFcns  called 

186 

187  %  Hint:  edit  controls  usually  have  a  white  background  on  Windows. 

188  %  See  ISPC  and  COMPUTER. 
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189  if  ispc 

190  setlhobject,  'BackgroundColor ' , 'white'); 

191  else 

192  setlhobject, 'BackgroundColor ', get (0, ' defaultUicontrolBackgroundColor ' ) ) ; 

193  end 

194 

195  % _  if 


196  %  -  Executes  on  mouse  press  over  axes  background. 

197  function  newname_ButtonDownFcn (hobject,  eventdata,  handles) 

198  %  hObject  handle  to  newname  (see  GCBO) 

199  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

200  %  handles  structure  with  handles  and  user  data  (see  GUI DATA) 

201 

202  % _  < 


203  %  -  Executes  on  button  press  in  pushbuttonl.  (Vector  Data  Manipulator  =  vdatamanipu 

lator ) 

204  function  pushbuttonl_Callback (hobject,  eventdata,  handles) 

205  %  hobject  handle  to  pushbuttonl  (see  GCBO) 

206  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

207  %  handles  structure  with  handles  and  user  data  (see  GUIDATA) 

208  global  filetoload  aflag  gflag  foldername  fitflag; 

209  if  (aflag==[])  aflag=0; 

210  end 

211  if  (gflag==[])  gflag=0; 

212  end 

213  if  (fitflag— []  )  fitflag=0; 

214  end 

215  %Determine  if  this  is  a  current  session  or  a  return  to  a  previous 
session  of  data  creation 

216  button  =  questdlg  (  '  Are  you  in  an  active  session  of  data  analysis? ',' Session  Type  Sele  if 
ction ' ) ; 

217  switch  button 

218  case  'Yes' 

219  nott='not'; 

220  [Valmanac, xnplot, ynplot , xinterp, yinterp, sample_name, BeamEnergy, probe_cur rent ]  if 
=vdatamanipulator ( filetoload,  aflag,  gflag,  foldername); 

221  almanac_save (Valmanac, foldername, xnplot , ynplot , xinterp, yinterp, gflag, nott , f il 
etoload) ; 

222  if  fitflag 

223  datatofit= [xinterp; yinterp] '; 

224  fbsLd (datatofit, foldername, sample_name, BeamEnergy, probe_current) ; 

225  end 

226  case  'No' 

227  prompt  =  {'Enter  the  folder  name  for  data  plot  storage Enter  name  of  superp  if 

osition  plot  to  amend'}; 

228  dlg_title  =  'Amended  Session  Input'; 

229  num  lines  =  1; 
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230  def  = 

231  answer  =  inputdlg (prompt , dlgtitle, num_lines , def ) ; 

232  answerstr=char (answer ) ; 

233  of=answer ( 1 ) ; 

234  oldf older=char (of ) 

235  superplotname=char (answerstr (2, :)); 

236  [ Valmanac, xnplot, ynplot, xinterp, yinterp, sample_name, BeamEnergy, probe_current] 
=vdatamanipulator ( f iletoload,  aflagl,  gflagl,  oldf older ) ; 

237  almanac_save (Valmanac, oldfolder, xnplot, ynplot, xinterp, yinterp, gflag, superplot 
name, f iletoload) ; 

238  if  fitflag 

239  datatofit= [xinterp; yinterp] '; 

240  fbsLd (datatof it , f oldername, sample_name, BeamEnergy, probe_current) ; 

241  end 

242  otherwise 

243  disp('To  find  the  superplot  and  folder  to  amend,  look  in  SEM\ .  Folder  data  se  i/ 
ts  are  named  by  a  lxLetter  designator  denoting  sample,  followed  by  the  beam  energy. 

The  superplot  is  in  this  folder.') 

244  end 

245 

246  % _ _  * 


247  %  -  Executes  on  button  press  in  pushbutton2  (ImageDataManipulator ) 

248  function  pushbutton2_Callback (hObject,  eventdata,  handles) 

24  9  %  hObject  handle  to  pushbutton2  (see  GCBO) 

250  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

251  %  handles  structure  with  handles  and  user  data  (see  GUIDATA) 

252  global  tiftoload  foldername  aflag  gflag  fitflag; 

253  if  (aflag==[])  aflag=0; 

254  end 

255  if  (gflag==[])  gflag=0; 

256  end 

257  if  (fitflag==[] )  fitflag=0; 

258  end 

259  %Determine  if  this  is  a  current  session  or  a  retail 
rn  to  a  previous  session  of  data  creation 

260  button  =  questdlg  (  '  Are  you  in  an  active  session  of  data  analysis? Session  Type  Sele  it 
ction ' , ' Cancel ' ) ; 

261  switch  button 

262  case  'Yes' 

263  nott='not'; 

264  [almanacR,  almanacC, xnplot, ynplot, xinterp, yinterp, sample  name, BeamEnergy, prob 
e_current] =imagedatamanipulator (tiftoload,  aflag,  gflag,  fitflag,  foldername); 

265  almanac_save  (almanacR,  foldername,  xnplot,  ynplot,  xinterp,  yinterp,  gflag,  nott,  tif  \t 
toload) ; 

266  flagc=0;  %Sets  the  superimposition  flag  to  'O'  so  column  data  * 

will  not  be  printed  on  superposition  plot 

2  67  almanac_save (almanacC, foldername, xnplot, ynplot,  xinterp, yinterp, f lagc,  nott,  tif 

toload) ; 
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268  if  fitflag 

269  datatofit=[xinterp;yinterp] '; 

270  fbsLd (datatof it , f oldername, sample_name, BeamEnergy, probe_current ) ; 

271  end 

272  case  'No' 

273  prompt  =  ('Enter  the  folder  name  for  data  plot  storage',  'Enter  name  of  superp  Z 

osition  plot  to  amend'}; 

274  dlg_title  =  'Amended  Session  Input'; 

275  num_lines=  1; 

276  def  =  (",''}; 

277  answer  =  inputdlg (prompt, dig  title, num_lines, def ) ; 

278  answerstr=char (answer) ; 

279  oldfolder= (answerstr (1, :)) ; 

280  superplotname=  (answerstr  (2,  : ) )  ; 

281  nott='not';  %Ensures  that  column  data  will  not  be 

282  [almanacR,  almanacC,  xnplot ,  ynplot ,  xinterp,  yinterp,  sample_name,  BeamEnergy,  prob  Z 
e_current] =imagedatamanipulator (tiftoload,  aflag,  gf lag, fitflag,  oldfolder) ; 

283  almanac_save (almanacR, f oldername, xnplot, ynplot, xinterp, yinterp, gf lag, superplo  Z 
tname, tiftoload) ; 

284  flagc=0;  %Sets  the  superimposition  flag  to  'O'  so  column  data  will/ 

not  be  printed  on  superposition  plot 

1 8 5  almanac_save (almanacC, foldername, xnplot, ynplot, xinterp, yinterp, f lagc, nott, tif  Z 

toload) ; 

286  if  fitflag 

287  datatofit= [xinterp; yinterp] '; 

288  fbsLd (datatof it, foldername, sample_name, BeamEnergy, probe_current) ; 

289  end 

290  otherwise 

291  displ'To  find  the  superplot  and  folder  to  amend,  look  in  SEM\ .  Folder  data  sets  ar Z 
e  named  by  a  lxLetter  designator  denoting  sample,  followed  by  the  beam  energy.  The  sZ 
uperplot  is  in  this  folder.') 

292 

293  end 

294 

295  % _  Z 


296  %  -  Executes  on  upcheck  of  checkboxl  =  This  is  for  'Keep  Almanac  for  mulitple  files Z 

297  function  checkboxl_Callback (hobject,  eventdata,  handles) 

298  %  hObject  handle  to  checkboxl  (see  GCBO) 

299  %  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 

300  %  handles  structure  with  handles  and  user  data  (see  GUIDATA) 

301  %  Hint:  get (hobject, 'Value' )  returns  toggle  state  of  checkboxl 

302  global  aflag; 

303  aflag  =  1; 

304 

305  % _  Z 


306  %  -  Executes  on  upcheck  of  checkbox2  for  super imposition  box. 


65 


C:\Documents  and  Settings\FMBradley\My  Documen. . . \newf irsttry.m 
November  16,  2005 


Page  9 
2:18:47  PM 


348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 


[legend_h, object_h, plot_h, text_strings]  =  legend(gca); 

[n, m] =size (text_strings ) ; 
counter=m/2; 

legendname=text_strings ; 
if  ishold 
yes=l; 

else  hold  on; 
end 

f irsttime=0; 

end 

fid  =  f open ( [ ' C : \Documents  and  Settings\FMBradley\My  Documents \ Physics \Data  AnalysisX  it 
' , folder, ' \Almanac . csv' ] , ' a ' ) ; 

fprintf (fid, '%14.12f  d  %d  %6.3f  %d  %6. 3f \n' , Dalmanac) ; 
f close (fid) ; 

if  firsttime==l  %print  superimposed  graphs  of  row  vector  data 

legendname={ } ; 

legendname=char (legendname) ; 

legendname=cellstr (legendname) ; 

counter=l; 

hold  on; 

grid  on; 

f irsttime=0; 

end 

if  flag  ==  1 

switch  counter 
case  {1} 

plot ( xnplot , ynplot , ' og ' ) ; 
plot (xinterp, yinterp, ' . g ' ) ; 

case  {2} 

plot (xnplot, ynplot, ’  ob' ) ; 
plot (xinterp, yinterp, ' .b ' ) ; 
case  {3} 

plot (xnplot, ynplot, ' oc ' ) ; 
plot (xinterp, yinterp,  ' . c  * ) ; 
case  {4} 

plot (xnplot, ynplot, 'ok' ) ; 
plot (xinterp, yinterp, ' . k ' ) ; 
case  {5} 

plot ( xnplot , ynplot , ' om ' ) ; 
plot (xinterp, yinterp, ' .m' ) ; 
case  {6} 

plot (xnplot, ynplot, ' oy ' ) ; 
plot (xinterp, yinterp, ' . y ' ) ; 
case  {7} 

plot (xnplot, ynplot, 'or') ; 
plot (xinterp, yinterp, ' . r ' ) ; 
case  {8} 

plot (xnplot, ynplot, ' sg ' ) ; 
plot (xinterp, yinterp, ' *g ' ) ; 
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function  [almanacR,  almanacC,  xnplot ,  ynplot,  xinterp,  yinterp,  sample_name,  BeamEnergy,  prob  if 
e_current] =imagedatamanipulator (tif toload, af lag, gf lag, f itf lag, folder) ; 

persistent  fileserial;  %Allows  var  fileserial  to  hold  value  if 

for  an  entire  session  of  multiple  calls  to  this  fn 

%if  loop  tests  if  variables  are  being  if 

saved  to  file, 

%and  increments  fileserial 
%to  help  differeniate 
%f igures 

f ileserial=f ileserial+1; 

%resets  all  figures  and  variables 
%if  ishold 
%  elf ; 

%  hold  off; 
lend 


%read  in  image  file 
10  =  imread  (tiftoload) ; 
10  =  double (10); 

10 ( : , 1) =0; 


Unitialize  variables 

Test=0; 

m=0 ; 


Test  for  location  of  the  spike  of  interest 


[ Y, I ] =max ( 10 ) ; 
Max) 

[K, L] =size  (I) ; 
MinVariance=300 ; 
.  m=l : (L-26) 
a=2  5 ; 


%Two  vectors:  Y=max  value  of  each  column,  I=Row#  of  each  columns  i f 
%K=row  size,  L=column  size 

festart  test  sequence  at  1  go  to  length  of  indice  matrix 

%test  for  remaining  testable  length  less  than  25  only  take  wh tf 


ats  left 

n=l:a  %Begin  creation  of  test  vector.  Populate  with  25  elements  of  I  (ri f 

ow  #  of  max  from  each  column) 

Test (n)=I ( (m+n) ) ; 


Variance=std  (Test )  ;  -'Calculate  variance  of  max  values'  row#s  to  select  desired  d  if 
ata  sample. 

if  Variance  <  MinVariance  lUpdate  minimum  variance 
MinVar iance=Variance ; 

seqstart=m;  %save  indice  start  point  in  vector  I  of  desired  data 

sample  (row  #s  of  desired  data  sample) 
seqstop=m+a; 

aeqlengr.h  ;  %  define  sequence  length  to  allow  for  seqlength  to  be  if 

used  as  indice  addition  term  for  creation  of  new  vector. 


end 
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end 

%Test  Row  data  for  least  variance  of  column  numbers  for  selected  column  data  sample 
I0T=I0  ' ; 

[ Y2 , 13] =max ( I0T) ;  %Two  vectors:  Y2=max  value  of  each  Row,  l3=Col#  of  each  row's  if 

Max) 

[KR,  LR] =size  (13) ; 

MinVarianceR=500; 
a=2  5 ; 

for  m3=2 : (LR-a-1)  %Loop  through  all  row  numbers 

for  n=l:a  %Begin  creation  of  test  vector.  Populate  with  25  elements  of  I  if 

2  (Column  #  of  max  from  each  row) 

TestR (n) =13 ( (m3+n) ) ; 

end 

VarianceR=std (TestR) ;  %Calculate  variance  of  maximum  values'  column  #s  to  if 

select  desired  data  sample. 

if  VarianceR  <  MinVarianceR  %Update  minimum  variance 
MinVarianceR=VarianceR; 

seqstartR=m3;  %save  indice  start  point  in  vector  12  of  desired  da  if 

ta  sample  (Column  #s  of  desired  data  sample) 
seqstopR=m3+a; 

seqlengthR=a;  %define  sequence  length  to  allow  for  seqlength  to  be 

used  as  indice  addition  term  for  creation  of  new  vector 
end 
end 

%Compare  MaxValue  vectors  and  choose  column  and  row  with  largest  same  maximum 
PeakPixelValue=0 ; 
for  stepR=l:25; 

for  stepC=l:25; 

if  Y (stepC+seqstart) ==Y2 (stepR+seqstartR) 

Peak=Y (stepC+seqstart) ; 
if  Peak>PeakPixelValue 
PeakPixelValue=Peak; 

MaxPixelColNum=seqstart+stepC; 

MaxPixelRowNum=seqstartR+stepR; 

end 

end 

end 

end 

%extract  row  and  column  data 
RowData=I0 (MaxPixelRowNum, : ) ; 

ColData=I0 ( : , MaxPixelColNum) ; 

%Create  Noise  vector  from  data  outside  of  spike 

%Calculate  variance  of  noise  and  through  out  sample  data  vectors 
for  z=l:50 

NoiseData (z) =RowData  ( z ) ; 
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85  end 

86 

87  RawVarOfNoise=std (NoiseData) ; 

88  RawMeanOfNoise=mean (NoiseData) ; 

89  for  zl=l: 50 

90  if  NoiseData (zl)  >=  RawMeanOfNoise+RawVarOfNoise 

91  NoiseData ( zl) =RawMeanOf Noise; 

92  end 

93  end 

94  MeanOfNoise=mean (NoiseData) ;  %Calculate  average  of  noise  within  variance 

95  VarOfNoise=std(NoiseData) ; 

96  %Normalize  row  vector  of  sample  data 

97  NormRowData= (RowData-MeanOfNoise) / (PeakPixelValue-MeanOfNoise) ; 

98 

99  %Normalize  ColData  vector 

100  for  z=l:50 

101  NoiseDataCol ( z ) =ColData ( z) ; 

102  end 

103 

104  RawVarOfNoiseCol=std(NoiseDataCol) ; 

105  RawMeanOfNoiseCol=mean (NoiseDataCol) ; 

106  for  zl=l:50 

107  if  NoiseDataCol (zl)  >=  RawMeanOfNoiseColt (2*RawVarOfNoiseCol) 

108  NoiseDataCol (zl) =RawMeanOfNoiseCol; 

109  end 

110  end 

111 

112  %Calculate  average  of  noise  within  variance 

113  MeanOfNoiseCol=mean (NoiseDataCol) ; 

114  VarOfNoiseCol=std (NoiseDataCol ) ; 

115  NormColData= (ColData-MeanOfNoiseCol) / (PeakPixelValue-MeanOfNoiseCol ) ;  i/ 


116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 


%add  spline  interpolation 

xr=  (MaxPixelColNum-49) : (MaxPixelColNum+50 ) ; 
xc=  (MaxPixelRowNum-49) : (MaxPixelRowNum+50) ; 
for  xl=l : 100 

yr(xl)  =  NormRowData (MaxPixelColNum-50+xl) ; 
yc (xl) =NormColData (MaxPixelRowNum-50+xl) ; 

end 

esc  =  spline (xc, [0  yc  0]); 
csr=spline (xr, [0  yr  0] ) ; 

xxr=linspace ( (MaxPixelColNum-50) , (MaxPixelColNum+50 ) , 1000 ) ; 
xxc  =  linspace ( (MaxPixelRowNum-50 ) , (MaxPixelRowNum+50 ) , 1000 ) ; 

%Extract  Halfmaxfullwidth  from  the  normalized  Row  data 
maxdif=l ;maxdif y=l ; 
for  xxx=l:500 
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133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 


dif=abs ( 0 . 2-ppval (csr , xxr (xxx) ) ) ; 
if  difcmaxdif 
maxdif=dif ; 
lhshalfmax=xxr (xxx) ; 

end 

end 

for  xxy=500 : 1000 

dify=abs (0 . 2-ppval (csr, xxr (xxy) ) ) ; 
if  dify<maxdify 
maxdify=dif y ; 
rhshalfmax=xxr (xxy) ; 

end 

end 

hold; 

Rhalfmaxfullwidth=rhshalfmax-lhshalfmax 

%Extract  the  halfmaxf ullwidth  from  the  normalized  column  data 
maxdi f  = . 2 ; maxdi f  y= . 2  ; 
for  xxx=l:500 

dif=abs (0 . 2-ppval (esc, xxc (xxx) ) ) ; 
if  dif<maxdif 
maxdif=dif ; 
lhshalfmax=xxc (xxx) ; 

end 

end 

for  xxy=500 : 1000 

dify=abs (0 . 2-ppval (esc, xxc (xxy) ) ) ; 
if  dify<maxdify 
maxdify=dify; 
rhshalfmax=xxc (xxy) ; 

end 

end 

hold; 

Chalfmaxfullwidth=rhshalfmax-lhshalfmax; 

%Assign  output  variables  to  almanac  vector 

prompt  =  ('Unique  sample  name ', 'Enter  probe  current  in  "3e-8"  notation: ',' Enter  Beam  / 
Voltage  in  format  "30"  kV:', 'Enter  Exposure  Time  in  "10.005"  (sec)  format:'}; 
dlg_title  =  'Almanac  File  Definition  Input'; 
num_lines=  1; 

def  =  {'S72','6e-10','25','1.0'}; 

answer  =  inputdlg (prompt, dlg_title, num_lines, def ) ; 

answerstr=char (answer) ; 

sample_name= (answerstr (1, : ) ) ; 

probe_current= (answerstr (2, : ) ) ; 

BeamEnergystring= (answerstr (3, : ) ) ; 

BeamEnergy=str2double (answerstr (3, :) ) ; 

ExpTime=str2double (answerstr (4,  : ) )  ; 

almanacR= [sample_name  str2double (probe_current )  0  BeamEnergy  Rhalfmaxfullwidth  PeakPi 
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xelValue  ExpTime] ; 

181  almanacC=  [sample_name  str2double (probe_current )  1  BeamEnergy  Chalfraaxfullwidth  PeakPi  if 
xelValue  ExpTime] ; 

182 

183 

184 

185  %  Page  1  graphs.  Overview  of  image  with  selected  max  values  in  colmax/row 

186  %  and  rowmax/col 

187  figure  (1) 

188  title ([' Overview  page  for:  '/folder,'  ', tif toload] ) ; 

189  subplot  (3, 1, 1) ,  imagesc  (abs  (10) )  ;  %Shows  reproduction  of  tiff  f  if 

ile  from  MicroCCD 

190  hold; 

191  plot (MaxPixelColNum, MaxPixelRowNum, ' xk ' ) ; 

192  title ([' MATLAB  reproduction  of  TIFF  (', tif toload, ')']) ; 

193  hold  off; 

194 

195  %  Plots  maximum  per  Row  vs  Column  number  where  that  maximum  falls  in  the  row. 

196  subplot  (3, 1, 2) ,  plot (max (10) ) ; 

197  hold; 

198  title (['Max  Pixel  value/Row  vs  Column  Number  of  MaxPixVal  (', tif toload, ')']) ; 

199  plot (MaxPixelColNum, (MeanOfNoise-VarOfNoise) ,  '  *y ' )  ; 

200  axis ( [ 0 , L, (MeanOfNoise-VarOfNoise ) , (max ( Y) +10 ) ] ) ; 

201  grid  on; 

202  hold; 

203 

204  %  Graph  Columns'  max  pixel  value  vs  row  Number  of  MaxpixVal 

205  subplot  (3, 1,  3) ,  plot  (max (I0T) )  ; 

206  hold; 

207  title (['Max  Pixel  value/Column  vs  Row  Number  of  MaxPixVal  (',  tif toload, ')']) ; 

208  plot (MaxPixelRowNum, (MeanOf NoiseCol-VarOfNoiseCol ) , ' *y ' ) ; 

209  axis([0,LR, (MeanOfNoiseCol-VarOf NoiseCol ) , (max (Y2) +10) ] )  ; 

210  grid  on; 

211  hold; 

212 

213  %Save  figure  1  to  file 

214  D=now; 

215  ddate=day  ( D )  ; 

216  [n,mmonth]  =month  (D)  ; 

217  dateserial=  [ '_' ,  num2str  (ddate)  , inmonth,  num2str  (file serial )  ]  ; 

218  filetosave= [ 'C: \Documents  and  Settings\FMBradley\My  Documents\Physics\Data  Analysis!' 

, folder , ' NOV' , sample_name, , num2str (BeamEnergy) , , probe_current, dateserial, ' . f ig ' 

]  ; 

219  saveas  (1,  filetosave)  ; 

220 

221  %Write  NormRowData  to  file 

222  Vdataf  ile=  [ '  C:  NDocuments  and  SettingsNFMBradleyNMy  DocumentsNPhysicsNData  Analysis!  ',  i/ 
folder,  ' !VRdata_' , sample_name, num2str (BeamEnergy) , '_' ,probe_current, dateserial, '.  if 
csv'  ]  ; 
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223  fid  =  fopen (Vdatafile, ' w' ) ; 

224  fprintf (fid, ' %11. 9f \n' , NormRowData) ; 

225  fclose(fid); 

226 

227  %Write  NormColData  to  file 

228  Vdataf ile= [ ' C : \Documents  and  Settings\FMBradley\My  Documents \Physics\ Data  Analysis'^', 
folder,  ' \VCdata_ ' , sample_name,  ' _ ' , num2str (BeamEnergy) ,  ,probe_current, dateserial, 
csv' ] ; 

229  fid  =  fopen (Vdataf ile, ' w ') ; 

230  fprintf (fid, ' %11 . 9f \n' , NormColData) ; 

231  fclose(fid); 

232 

233  %Plot  Raw  row  and  column  data  sets 

234  figure  (2 )  ; 

235  title (['Peak  plots  for:  ', folder tiftoload] ) ; 

236  %Raw  Row  and  Column  data 

237  subplot (2, 2, 1) ,  plot (RowData, ' o ' ) ; 

238  title (['Raw  Data  for  Peak  Pixel  Row  (', tiftoload, ')']) ; 

239  axis ( [ (MaxPixelColNum-50) , (MaxPixelColNum+50 ) , 0, (PeakPixelValue+10) ] ) ; 

24  0  grid  on; 

241  subplot (2, 2, 2 ) ,  plot (ColData, 'o ' ) ; 

242  title  (['Col  Data  for  Peak  Pixel  Column  (', tiftoload, ')']) ; 

243  axis ( [ (MaxPixelRowNum-50) , (MaxPixelRowNum+50 ) , 0, (PeakPixelValue+20) ] ) ; 

24  4  grid  on; 

245 

246  %Print  Normalized  Row  Data 

247  R=0;xRlimit=0;xRstep=0; 

248  [R,xRlimit]  =size (NormRowData) ; 

24  9  for  xRstep=l:xRlimit; 

250  XRnorm(xRstep) =xRstep-MaxPixelColNum; 

251  end 

252  subplot  (2, 2,  3) ,  plot  (XRnorm,  NormRowData,  1  o 1 )  ; 

253  hold; 

254  plot ( (xxr-MaxPixelColNum) , ppval (csr, xxr) 

255  title ([ 'Normalized  Row  Data  for  Peak  Pixel  Row  (', tiftoload, ')')) ; 

256  axis(  [(-50) ,  <50)  ,-.2, 1.1] )  ; 

257  grid  on; 

258  text (10, 0.5, [ ' FWHM^ ' , num2str (Rhalfmaxfullwidth) ] ) ; 

259 

260 
261 
262 

263  %Print  Column  Data 

264  [xlimit,  C]  =size  (NormColData)  ; 

265  for  xstep=l : xlimit; 

266  XCnorm (xstep) =xstep-MaxPixelRowNum; 

267  end 

268  subplot  (2, 2,  4) ,  plot  (XCnorm,  NormColData,  '  o ')  ; 

269  hold; 
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270  plot ( (xxc-MaxPixelRowNum) ,ppval (esc, xxc) 

271  title ([' Normalized  Col  Data  for  Peak  Pixel  Column  ( ' , tif toload, ' ) ' ] ) ; 

272  axis ([(-50), (50), -.2,1.1]); 

273  grid  on; 

274  text (10,0.5, [ ' FWHM= 1 , num2str (Chalfmaxfullwidth) ] ) ; 

275 

276  %Save  figure  2  to  file 

277  f iletosave2= [ ' C : \Documents  and  Settings\FMBradley\My  Documents\Physics\Data  AnalysisN 
', folder, ' \RC_' , sample_name, 1 _ ' , num2str (BeamEnergy) , ' _ 1 , probe_current,  dateserial, 

fig' ] ; 

278  saveas(2,  filetosave2) ; 

279 

280  %Test  for  desire  for  superimposed  graph  page  and  set  export  variables  for 

281  %consolidation  in  newfirsttry 

282  if  (or  (gflag,  fitf lag)  ) 

283  xnplot=XRnorm; 

284  ynplot=NormRowData; 

285  xinterp=xxr-MaxPixelColNum; 

286  yinterp=ppval  (csr, xxr )  ; 

287  else 

288  xnplot=0; 

289  ynplot==0; 

290  xinterp=0; 

291  yinterp=0; 

2  92  end 
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1  % - Created  by  FM  Bradley  to  perform  a  least  squares  fit  for  data  to  the 

2  %2-D  model  for  diffusion  of  minority  carriers  in  a  luminscent 

3  %semiconductor . 

4 

5  %Revision  Log:  13  September  -  Time  remaining  window  added  by  FMBradley 

6 

7  %Takes  a  2-D  vector  data  in  (xstepX2)  size,  normalizes,  cuts  to  1%  of  max, 

8  %performs  a  least  squares  fit  to  a  2-D  model  of  carrier  distro  given  a 

9  %guassian  distributed  generation  area.  Output  is  the  fitted  model  vector, 

10  %the  difference  between  model  and  experimental  data,  plots  of  the  parameter 

11  %space,  both  curves,  the  selected  of  optimal  Mu,  and  n=2sigma  radius  of 

12  %generation  area  that  fits  the  data. 

13 

14  function  [FitteaCurve, residue, Paramspace, Lopt, nopt ] =fbsLd (datain, folder, sample_name, B 
eamEnergy, probe_current ) 

15 

16  %CreateStruct . WindowStyle= ' replace ' ; 

17  %  CreateStruct .  Interpreter '  tex ' ; 

18  %  h=msgbox (' Calculating . ','Progress  Monitor ', CreateStruct ) ; 

-  '  .  -  :  %Significant  radius  of  generation  volume  distribution (cm)  95%  of  cur re  i/ 

nt  ring?? 

20  V=0;  ^Applied  bias  (Volts) 

21  %tau=4 . 32e-9;  %input  variable.  Lifetime  of  carrier  in  question  (s) 

22  icd=1.0;  %inter-contact  distance  (mm)  for  calculation  of  electric  field  assoc  if 

iated  with  V 

23  E=V/icd*10;  Electric  field  (V/cm) 

24  step=0 . 04e-4 ;  %step  size  in  (cm)  designed  to  fit  the  incr 

ement  size  of  our  experimental  data  w/interpolation 

25  [R3, 13] =size (datain) ;  %Id  dimension  of  datain  array 

2  6  datain {:, 1) =0 . 4e-4 . *datain (:,  1 } ;  %Use  this  line  if  the  data  is  coming  from 

actual  experimentally  gathered  and  spline  interpreted  data  from  the  MATLAB  GUI. 

27  %Necessary  because  spline  data  is  already  if 

in  microns  and  we  must  convert  all  numbers  here  to  cm  if 

28  Lmin=SlopeL2 (datain) -. 3e-4 ;  %Conducts  In  slope  analysis  method  for  t if 

he  diffusion  length  assessment 

29  Lstep=.le-4; 

30  stop=0; 

31  while  (stop<3) 

32  switch  stop 

33  case  0 

34  Lstep=.2e-4; 

35  nstep=.2e-4; 

36  Lmax=Lmin+ . 4e-4 ; 

37  nmax=nmin+ . 4e-4  ; 

38  case  1 

39  Lstep=. le-4; 

40  nstep=.le-4; 

41  Lmax=Lmin+ . 2e-4 ; 
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42 

nmax=nmin+ . 2e-4 ; 

43 

case  2 

44 

Lstep=. le-4; 

45 

nstep=. le-4; 

46 

nmin=mnin- . 5e-4  ; 

47 

Lmin=Lmin- . 5e-4 ; 

48 

Lmax=Lmin+le-4 ; 

49 

nma  *nmin+ le-  4  ; 

50 

end 

51 

52 

a=0 

;  %step  variable  for 

indexing  x 

53 

b=0 

■  tstep  variable  for 

indexing  y 

54 

MinSum=10000; 

55 

istep=0; 

56 

Paramspace=0 ; 

57 

for 

L=Lmin : Lstep : Lmax 

58 

t0=cputime; 

59 

n=nmin; 

%1-d  spot  dimension  of  a  squa 

re  generation  area  (cm) 

60 

istep=istep+l; 

%parameter  space  indexing  var  4 

iable  for  mu 

61 

or  n 

jstep=0; 

%param  space  index  variable  f  t/ 

62 

S= (E*LA2) / . 025 ; 

%Field  effect  (cm) 

63 

mtproduct=LA2/ .025; 

%mu  tau  product  in  (cm/'2/V) 

64 

for  n=nmin:nstep:nmax 

65 

G=2/n; 

%Gaussian  2sigma  radius  f 

or  95%  generation  distro  (unitless) 

66 

j  step= j  step+1 ; 

67 

xmin=datain (1,1) ; 

^convert  datain  from  pixels  to  cm  de 

pendent 

upon  resolution  size  of  SEM/Optical  scope  . 4microns/pixel 

68 

xmax=max ( datain ( : ,  1 )  )  ; 

69 

%  (cm) 

x=xmin; 

* 

70 

%  (cm) 

ymin=0; 

</ 

71 

%  (cm) 

ymax=0 ; 

* 

72 

y=ymin; 

73 

b=0; 

74 

for  yl=ymin: step: ymax 

75 

a=0; 

76 

b=b+l; 

77 

for  xl=xmin: steprxmax 

78 

lntegrand=0; 

79 

a=a+l; 

80 

Integrand=dblquad ( Qlntgrnd,  (xl-n) ,  (xl+n) ,  (yl-n) ,  (yl+n) ,  step,  [ ] , xl, yl, 

S,L,G)  ; 

81 

if  (isnan (Integrand) ) 

Integrand=OldIntegrand+ . 01*01dlntegrand; 
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82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 


end 

01dlntegrand=lntegrand; 

Int (a,b, l)=xl; 

Int (a, b, 2 ) =yl; 

Int  (a,  b,  3)  =  (G/ (2‘piA2+LA2) ) ‘Integrand; 
end  %end  x  loop 

end  %end  y  loop 
%Least  square  computation 

[M2, I2]=max (Int ( : , 1, 3) ) ; 

NormInt= [Int ( : , 1, 1 )  (Int ( : , 1, 3) . /M2) ] ; 
%LogNormInt=log (Normlnt {: ,2) ) ; 

%LogDataIn=log (normdatain ( : , 2) ) ; 

LS=0  ; 


for  loopc=l:R3;  %1 :  (  (  (dimof data-1 ) /2 ) -105 )  ;  W 

%start  loop  at  max  +13  and  count  to  end  testing  for  least  squares  diff 
%LS (loopc) = (LogNormlnt ( I2+100+loopc) -LogDataln (13+100+loopc) ) A2 
%;  Used  to  calculate  log  difference  squares  between  two  lines 
%of  data 

if  isnan (Normlnt (loopc, 2) ) 

Normlnt (loopc, 2) =1; 

end 

if  isnan (datain (loopc, 2 ) ) 
datain (loopc, 2 ) =1; 

end 

if  datain (loopc, 2) <. 7 

he  top  30%  of  the  curve  that  does  not  fit  our  assumptions 

LS (loopc)  = (Normlnt (loopc, 2 ) -datain (loopc,  2) ) A2; 
end 

el  fit  data 
end 

quares  sum  loop 

LoopSum=sum (LS) ; 

Paramspace (istep, jstep, 1) = (istep-1) ‘Lstep+Lmin; 

Paramspace (istep, jstep, 2) = ( jstep-1) ‘nstep+nmin; 

Paramspace (istep, jstep, 3) =LoopSum; 
if  LoopSum<MinSum 
MinSum=LoopSum; 
residue=MinSum; 

ModelFit=NormInt; 

%Muopt=mu; 
nopt=n; 

Lopt=L; 
imin=istep; 
jmin=jstep; 

end 

figure (7);  %Plot  the  search  space  each  round  of  calc 

ulations  complete 

nplot=n*le4; 


%Eliminates  t 

%ends  non  modrf 
lends  Least  s  i/ 
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126 

127 

128 

129 

130 

131 

132 

133 

134 

135 


136 
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139 

140 
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143 
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147 

148 
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150 
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160 
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172 


Lplot=L*le4 ; 

if  (stop==l ) plot (Lplot , nplot, ' xr ' ) ; 

else  plot (Lplot, nplot, ' xk ') ; 

end 

if  not(ishold)  hold; 
end 

end  %end  of  n  for  loop 

%  CreateStruct . WindowStyle= ' replace ' ; 

%  CreateStruct .  Interpreter '  tex  * ; 

%  h=msgbox ([' Total  time  remaining=  ' , num2str ( ( (cputime-tO ) /60 ) * (Lmax-L) /Lstep) , '  m 
inutes.  Step  ' , int2str (istep* jstep) , '  of  ' , num2str (({ (Lmax-Lmin) /Lstep) +1 )*(( (nmax-nm 
in) /nstep) +1) ) , '  Total  Steps '],' Progress  Monitor ', CreateStruct) ; 
end  %end  of  L  for  loop 

%Decision  tree  for  creating  direction  of  propogation  of  parameter  space 
%build  up. 
switch  imin 
case  1 

Lmi n=Lmi n - L s t ep ; 
if  jmin==l 

nmi n=nmi n - n s t ep ; 
elseif  jmin==2 
nmin=nmin; 

else  nmin=nmin+nstep; 
end 
case  2 

Lmin=Lmin; 
if  jmin==l 

nmin=nmin-nstep; 
elseif  jmin==3 

nmin=nmin+nstep; 

else 

stop=stop+l; 

plot ( (Lmin*le4) , (nmin*le4) , '  xb' ) ; 

end 
case  3 

Lmin=Lmin+Lstep; 
if  jmin==l 

nmin=nmin-nstep; 
elseif  jmin==2 
nmin=nmin; 

else  nmin=nmin+nstep; 
end 

otherwise 

break; 

end  %end  for  the  switch  loop 

end  %end  for  the  while  loop 

title (' Search  Space');  Label  Figure  7 

xlabel ( ' Ld  ( \mum) ' ) ; 
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173  ylabel (' 2\sigma  (\mum)  ' )  ; 

174  hold;  %Hold  off  for  figure  7 

175  FittedCurve=ModelFit ; 

176  %Write  Fitted  Curve  data  to  file 

177  Vdataf ile= [ ' C : \Documents  and  Settings\FMBradley\My  Documents\Physics\Data  Analysis!',  it 
folder, ' \ModelFit ' , sample_name,  ' Ld_ 1 , Lopt,  ' _ ' , num2str (BeamEnergy) ,  ' kV_ ' , probe_current  it 
, ' .csv' ] 

178  fid  =  fopen (Vdatafile, 'w' ) ; 

179  fprintf  (fid,  '  %11 . 9f  \n ' ,  ModelFit )  ; 

180  fclose(fid); 

181  Residue=MinSum 

182  nopt 

183  Lopt 

184 

185 

186  %plot  the  resulting  fitted  curves 

187  figure  ( 8 )  ; 

188  plot (ModelFit ( : , 1 ) , ModelFit ( : , 2 ) ,  ' .b ' , datain ( : , 1) , datain ( : , 2 ) ,  ' -r  ' )  ; 

189  hold  on; 

190  Xmax=xmax  ; 

191  Xmin=xmin; 

1 92  Yrain=0 ; 

193  Ymax=1.0; 

194  Axis ( [Xmin  Xmax  Ymin  Ymax]); 

195  title ( 'Normalized  Experimental  vs.  Model  Fit'); 

196  xlabel (' Radial  distance  from  beam  center  (cm)'); 

197  ylabel ( 'Normalized  Intensity'); 

198  % ’ \mum: \mu= ' , num2str (Muopt) , ' : \tau= ' , num2str ( tau*le9)  , 

199  legendnamel=  [ 'Model  Fit :  (n= ',  num2str  ( nopt  *le4 ),'  \mum:  Ld= ' , num2str (Lopt *le4 ) , ' \mum)  R  it 
esidue=' , num2str (residue) ] ; 

200  legendname2= [ ' Data :  ' , sample_name,  ' : PC= ' , num2str (probe_current) , 'Amps ' ] ; 

201  legend(legendnamel, legendname2) ; 

202  hold  off; 

203 

204 

205  %Plot  the  parameter  space  of  the  resulting  IMicronXIMicron  parameter  range  for  both  2  it 
sigma  and  Ld 

206  figure  (9); 

207  contour (Paramspace ( : , : , 2) , Paramspace ( : , : , 1) , Paramspace ( : , : , 3) , 100) ; 

208  hold  on; 

209  plot  (nopt,  Lopt,  '  *y' ) 

210  Xmax=nmax  ; 

211  Xmin=nmin; 

212  Ymin=Lmin; 

213  Ymax=Lmax ; 

214  Axis ( [Xmin  Xmax  Ymin  Ymax] ) ; 

215  title  (['  Parameter  Space:  ',  sample_name,  num2str  (BeamEnergy)  ,'  kV:  ',  num2str  (probe_cu  it 
rrent) , ' A ' ] ) ; 

216  ylabel('Ld  step  multiples  [cm/'2/Vs]  ' )  ; 
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217  xlabel ( ' 2\sigma  steps  of  generation  area  [cm] ' | ; 

218  colorbar  ( '  vert ' )  ; 

219  hold  off; 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
■7 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 


function  [almanacV, xnplot , ynplot , xinterp, yinterp, sample_name, BeamEnergy, probe_current 
] =vdatamanipulator (f iletoload,  af lagl,  gflagl,  foldername)  ; 


persistent  fileserial;  %Allows  var  fileserial  to  hold  value 

for  an  entire  session  of  multiple  calls  to  this  fn 

%if  loop  tests  if  variables  are  being 


saved  to  file, 


%and  increments  fileserial 
%to  help  differeniate 
%f igures 

fileserial=fileserial+l; 


%resets  all  figures  and  variables 
%if  ishold 
%  elf ; 

%  hold  off; 

%end 

%Row  data  input  manipulation  loop 
%filename=input ( 'Enter  row  data  file  name: ' , 's ' ) ; 
rowdatainput=load ( f iletoload) ; 

RowData=rowdatainput  ( : , 2)  ; 

%Data  preparation  loop 

[PeakPixelValue, MaxPixelColNum] =max (RowData) ; 

%Create  Noise  vector  from  data  outside  of  spike 

%Calculate  variance  of  noise  and  through  out  sample  data  vectors 
for  2=1:200 

NoiseData ( z ) =RowData ( z )  ; 

end 


RawVarOfNoise=std (NoiseData) ; 

RawMeanOfNoise=mean (NoiseData) ; 
for  zl=l : 200 

if  NoiseData ( zl )  >=  RawMeanOfNoise+RawVarOfNoise 
NoiseData ( zl) =RawMeanOf Noise  ; 

end 

end 

..  fM  ise=mean(NoiseDat  ;  % Calculate  average  of  noise  within  variance 

VarOfNoise=std (NoiseData) ; 

%Normalize  row  vector  of  sample  data 
NormRowData= (RowData-MeanOfNoise) / (PeakPixelValue-MeanOfNoise) ; 


%add  spline  interpolation 

xr=  (MaxPixelColNum-49) : (MaxPixelColNum+50 ) ; 
for  xl=l : 100 

yr(xl)  =  NormRowData (MaxPixelColNum-50+xl ) ; 

end 
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47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
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62 
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64 

65 

66 

67 

68 
69 
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71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 


csr=spline (xr, [0  yr  0]  )  ; 

xxr=linspace ( (MaxPixelColNum-50) , (MaxPixelColNum+50 ) , 1000 ) ; 

%Extract  Halfmaxfullwidth  from  the  normalized  Row  data 

maxdif=l;maxdify=l; 

for  xxx=l:500 

dif=abs (0 . 2-ppval (csr , xxr (xxx) ) ) ; 
if  difcmaxdif 
maxdif=dif ; 
lhshalfmax=xxr (xxx) ; 

end 

end 

for  xxy=500 : 1000 

dify=abs (0 . 2-ppval (csr, xxr (xxy) ) ) ; 
if  difycmaxdify 
maxdify=dify; 
rhshalfmax=xxr (xxy) ; 

end 

end 

Rhalfmaxfullwidth=rhshalfmax-lhshalfmax; 

%Assign  output  variables  to  almanac  vector 

prompt  =  ('Unique  sample  name', 'Enter  probe  current  in  "3e-8"  notation: ', 'Enter  Beam  i/ 

Voltage  in  format  "30”  kV:', 'Enter  Exposure  Time  in  "10.005"  (sec)  format Enter  R 

or  C  for  Row  or  Column  Data'); 

dlg_title  =  'Almanac  File  Definition  Input'; 

num_lines=  1; 

def  =  {'S72','6e-10\ '25\ '1.0','R'}; 
answer  =  inputdlg (prompt, dlg_title, num_lines, def ) ; 
answerstr=char (answer) ; 
sample_name= (answerstr (1,  :)); 
probe_current=(answerstr (2,  :) ) ; 

BeamEnergystring= (answerstr (3,  : ) ) ; 

BeamEnergy=str2double (answerstr (3,  : ) ) ; 

ExpTime=str2double (answerstr ( 4,  : ) ) ; 
switch  answerstr (5) 
case  'R' 

vector_type= 'O'; 
case  'C' 

vector_type=' 1 ' ; 
otherwise 

dialgl'Your  test  did  not  work'); 

end 

almanacV= [str2double (probe_current)  vector_type  BeamEnergy  Rhalfmaxf ullwidth  PeakPixe 4 
lvalue  ExpTime] ; 

%Analysis  output  loop 

%Plot  Raw  row  and  column  data  sets 

figure ( 1 ) ; 
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93  %Plot  overview  graph  of  row 

94  subplot  (1, 2, 1) ,  plot (RowData, 'o' ) ; 

95  title (['Raw  Data  for  Peak  Pixel  Vector  (', foldername, f iletoload, ')'])  ; 

96  axis ( [ (MaxPixelColNum-50) , (MaxPixelColNum+50 ) , 0, ( PeakPixelValue+10 ) ] ) ; 

97  grid  on; 

98 

99  %Print  Row  Data 

100  [xRlimit, R] =size (NormRowData) ; 

101  for  xRstep=l : xRlimit ; 

102  XRnorm (xRstep) =xRstep-MaxPixelColNum; 

103  end 

104  subplot (1, 2, 2) ,  plot (XRnorm, NormRowData,  '  o 1 )  ; 

105  hold; 

106  plot ( (xxr-MaxPixelColNum) , ppval (csr ,  xxr )  ,'.'); 

107  title ([ 1  Normalized  Data  for  Peak  Pixel  Vector  (',  f iletoload, ')']) ; 

108  axis([(-50),  (50), -.2,1.1]); 

109  grid  on; 

110  text (10,0.5, [ ' FWHM^ ' , num2str (Rhalfmaxfullwidth) ] ) ; 

111  hold; 

112 

113  %Save  figure  1  to  file 

114  D=now; 

115  ddate=day  (D)  ; 

116  [n,mmonth]  =month  (D)  ; 

117  dateserial=  [ '  ,  num2str  (ddate)  , inmonth,  num2str  (fileserial)  ]  ; 

118  filetosave=  [ 'C:  NDocuments  and  Settmgs\FMBradley\My  Documents\Physics\Data  AnalysisVtf 

,  folder  name,  '  \  ' ,  answer  str  (4 )  ,  num2str  (BeamEnergy)  ,  ,  probecurrent,  dateserial,  '  .  f  ig '  ] 

119  saveas  (1,  filetosave)  ; 

120  figure  (2); 

121 

122  %Set  export  variables  for 

123  %consolidation  in  newfirsttry 

124  xnplot=XRnorm; 

125  ynplot=NormRowData; 

126  xinterp=xxr-MaxPixelColNum; 

127  yinterp=ppval (csr , xxr) ; 
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APPENDIX  E.  SLOPE  ANALYSIS  ALGORITHM  (SLOPEL2 . M) 


C:\Documents  and  Settings\FMBradley\My  Documents\P. . . \SlopeL2 .m  Page  1 

1  function  [L] =SlopeL (datavector) 

2  % 

3  % 

4  % - Automated  slope  determination  algorithym 

5  %to  gain  initial  indication  of  diffusion  length  of  a  material. 

6  %  Original  protocol  determined  by  Dave  Luber.  Automated  by  Mitch  Bradley 

7  %  3  October,  2005 

8  j=l; 

9  y=0; 

10  i-1; 

11  [V, I] =max (datavector ( : , 2) ) ; 

12  EOV=size (datavector (:,  1) ) ; 

13  for  i=l:I 

14  y=datavector (i,  2)  ; 

15  if(yc.l)  %These  limits  should  be  adjusted  for  noise  level  of  the  sai/ 

mple  image  (200M  above  lower  limit) 

16  if  (y>.01)  %This  should  be  adjusted  for  lower  limit  of  noise  leve  if 

1  for  sample  image 

17  yfit ( j ) =log (y)  ; 

18  xfit ( j ) =datavector (i,  1)  ; 

19  j=j+lf 

20  end 

21  %i=i+l; 

22  end 

23  end 

24  p=polyfit (xfit, yfit, 1) ;  iconduct  polylinomial  fit  for  poly  of  degree  if 

1  (linear)  for  data  in  xfit  and  yfit 

25  L=l/p  (1) 

26  b=p ( 2 ) 

27  fitdata= [xfit ' ,  yfit1]; 

28 

2  9  % _ Loop  for  the  right  side  of  the  distribution _ if 


30  j-1; 

31  for  i=I : EOV 

32  yr=datavector (i,  2 )  ; 

33  if(yrc.l)  %These  limits  should  be  adjusted  for  noise  level  of  the  s  if 

ample  image  (200M  above  lower  limit) 

34  if  (yr>.01)  %This  should  be  adjusted  for  lower  limit  of  noise  lev  if 

el  for  sample  image 

35  yfitr  (j)=log(yr)  ; 

36  xfitr (j ) =datavector (i, 1) ; 

37  j-j+1; 

38  end 

39  %i=i+l; 

4  0  end 

4 1  end 

42  pr=polyf it (xfitr , yfitr,  1) ;  %conduct  polylinomial  fit  for  poly  of  degr  if 

ee  1  (linear)  for  data  in  xfit  and  yfit 
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4  3  Lr=-l/pr(l) 

44  br=pr  (2) 

45  fitdataR= [xfitr ' ,  yfitr ’ ] ; 

46 

47  % - Plot  In (Intensity)  vs  X  pos  for  all  data  left  of  spot - 

4  8  figure (4); 

4  9  plot (datavector ( : , 1) , log (datavector ( : ,  2 )  )  )  ; 

50  %axis ( [ (datavector (1,1)), (0) , (log (min (datavector ( : , 2) ) ) ) , (0) ] ) ; 

51  title ( 'LN( Intensity)  of  Left  Side  Distribution'); 

52  xlabel (' Radial  distance  from  beam  center  (cm)'); 

53  ylabel ( 'LN (Normalized  Intensity)'); 

54 

55  % - Calculate  error  reports  and  plot  line  slopes - 

56  figure  (5)  ; 

57  Lplus=L+ ( . 4e-4) / (max (yf it ) -yf it (1) ) ;  %The  Plus  factor  should  be  .4e-4  for 

real  data  and  . 04e-4  for  tfbs  control  data 

58  Lminus=L- ( . 4e-4 ) / (max (yf it ) -yf it ( 1 ) ) ;  %The  minus  factor  should  be  . 4e-4  for 

real  data  and  .04e-4  for  tfbs  control  data 

59  plot  (xf  it,  yf  it,  '  ob' )  ; 

60  hold; 

61  plot (xf it, polyval (p, xf it ) , ' -r ' ) ; 

62  yplus=(l/Lplus) *xfit+b; 

63  yminus= (1/Lminus ) *xf it+b; 

64  plot (xfit, yplus, ' -g ') ; 

65  plot (xfit, yminus, ' -k' ) ; 

66  % _ Right  side _ 

67  Lplusr=Lr+ (. 4e-4) / (max (yfitr ) -min (yfitr) ) ;  %The  Plus  factor  should  be  .4e-tf 

4  for  real  data  and  .04e-4  for  tfbs  control  data 

68  Lminusr=Lr- (.  4e-4)  /  (max  (yfitr ) -min  (yfitr ))  ;  %The  minus  factor  should  be  .  4e  \£ 

-4  for  real  data  and  .04e-4  for  tfbs  control  data 

69  plot (xfitr, yfitr, ' ob' ) ; 

70  plot  (xfitr,  polyval  (pr,  xfitr ),'  -  r  ')  ; 

71  yplusr= (-1/Lplusr ) *xf itr+br ; 

72  yminusr= (-1/Lminusr) *xfitr+br; 

73  plot (xfitr , yplusr, ' -g ') ; 

74  plot (xfitr , yminusr,  ' -k ') ; 

75 

76  legendnamel= [ ' Data  Points']; 

77  legendname2= [ ' L  actual  Ld= ' ,  num2str (L*le4 ) ,  ' \mum ' ] ; 

78  legendname3=[ 'L  Ld(  +  )  ■= ' ,  num2str  (Lplus*le4) ,  '  \mum'  ]  ; 

79  legendname4= [ ' L  Ld (-)  =  ',  num2str (Lminus*le4 ),' \mum' ] ; 

8  0  legendname5= [ ' Data  Points ' ] ; 

81  legendname6= [ ' R  actual  Ld=' , num2str (Lr*le4) , ' \mum' ] ; 

82  legendname7= [ ' R  Ld (+)=', num2str (Lplusr*le4 ),' \mum' ] ; 

83  legendname8=[ 'R  Ld (-)=', num2str (Lminusr*le4) ,' \mum' ] ; 

84 

85 

8  6  legend (legendnamel, legendname2, legendname3, legendname4, legendname5, legendname6, legend 
name7, legendname8) ; 
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87  title ('Slope  Analysis  Estimate  of  Ld' ) ; 

88  xlabel (' Radial  distance  from  beam  center  (cm)'); 
8  9  ylabel  ( 1 LN  (Normalized  Intensity)'); 

90  hold; 

91 
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APPENDIX  F.  NUMERICAL  INTEGRATION  SOLUTION  FOR 
MINORITY  CARRIER  DISTRIBUTION  ( INTEGRAND . M) 
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1  function  I=Intgrnd(chi,  ada,  xl,  yl, S, L, G) 

2  I=exp (- (GA2) .*( (xl-chi) . A2+(yl-ada) . A2 ) + (S . *chi/ ( 2*L . A2 ) ) ) . *besselk(0,  (sqrt ( (S. A2+4*L 
. A2) .* (chi. A2+ada. A2) ) / (2*L. A2) ) )  ; 

3  %I=exp  (S .  *chi/  (2*L.  A2 )  )  .  *besselk  (0,  (sqrt  (S.  A2+4*L.  A2)  .  *  (chi  .  A2+ada.  A2)  .A(.5)/(2*L.A2)  <£ 

)  ) ! 

4 
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APPENDIX  G.  PHOTON  RECYCLING  PERTURBATION 

REFERENCE  EQUATION  SHEET 


Calculated  Properties 


Measured  Pritfisr 

Na=i 


Referenced  Properties 


01  :=1 .21  Rt=,001 
ffi:=1.2l  R2=.00l 
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SAMPLE  D6 
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